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Section 1 


Two Dimensional Subsonic Flow 


of Compressible Fluids 


Summary 


The basic concept of the present peper is 
to use the tangent line to the adiebetic pressure- 
-volume curve es an approximation to the curve 
itself. First, the general characteristics of 
such ^X fluid are shown. Then in Section I,a 
theory is developed Which in — ic 
similar to thet af DemtchXgenko sn? Buserenn 
but is more general and can be applied to flow 
with velocity approach) thet of soe ]] The 
theory is then applied to celerlate the flow 
over elliptic cylinders. In Section II the 
work of h. Retemen is applied to thie epprovimatejd 
ediabatic fluid and &wo-eguaiione-ame—piven— which 


speed charatonistice of sinfoile. 
th e e necati. LA 


TWO-DIMENSIONAL SUBSONIC FLOW 


OF COMPRESSIBLE FLUIDS 
ty 


Heue-shen Tsien 
California Institute of Technology 


Introduction 


Assuming that the pressure is a single=valued function of 
density only, the equations of two-dimensional irrotational motion 
of compressible fluids can be reduced to a single non-linear equation 
of the velocity potential. In the supersonic oase, that is, in 
the case when te flow velocity is everywhere greater than that of 
local sound velocity, the problem is solved by Meyer Prandtl and 
Busemann using the method of characteristics. The essential difficulty 
of this problem lies in the subsonic case, that is, in the case when 
the flow velocity is everywhere smaller than but near te the local 
sound velocity, because then the method of characteristios cannot 
be used, Glauert & Prandtl (Ref. 1) treated the case when the 


disturbance to the parallel rectilinear flow due to presence of a 


solid body is small, They were then able to linearize the differential 
equation for the velocity potential and obtained an equation very 
similar to eme for Né incompressible fluids. But there are usually 
stagnation points either in the surface of the body or in the field 
of flow, where the disturbance is no longer small, Hence, it is 
doubtful whether the linear theory oan be applied to the flow near a 
stagnation point, On the sane ground, the theory breaks dom in the 
ease of bodies whose dimension across the stream is not small compared 
with the dimension parallel to the stream. 

— . SRR, Janzen 
and Rayleigh developed the method of successive approximations. This 


«nte 
method A into a more convenient form 
by L. Poggi G and P.A. Walther {Met}, Recently C. Kaplan 
(Ref. 8) treated the case of flow over Joukowsky airfoils and elliptic 
cylinders, using Poggi's method, However, the method is rather tedious 
and the convergent very slow if the local velocity of sound is approached, 
"Luck" Molenbroek (= and Tschapligin =) suggested 
the use of the magnitude of velocity 4/ anà inclination Z of 
vend chosen axis as independent variables, and were thus able 
to reduce the equation of velocity potential to a linear equation. 
This equation was solved by Tschapligin, (IF) emirenmtty-put 
PPP ⁵² 8 ieasen—md-M Sinmmer (amt). 
The solution is essentially a series each term of which is a product 
of a hypergeometrio function of 4 end a trigonometric function 
of Ø . The min difficulty in practical application of this 
solution is to obtain a proper set of boundary conditions in the 


plane of independent variables wp and put the solution in a 


closed form. l 

Techapligin e shows that a great simplification of 
the equation in the hodograph plane results if the ratio of the specific 
heats of the gas is equal to -l — , Since all real gases have their 
ratio of specific heats between l and 2, the value -l seems without 
practical significance. It was Demtdohenko (Ref. 3) ma Duseman 
(Ret, 4) who made Acar, the moaning of this specific value of -1. 
They found that this really means to take the tengent of, pressure- 
volume curve as an approximaté to the curve itself. However, they 
limit themselves to use the tangent at the state of rt the gas — 


Thus their theory oan only be applied to flow with velocities up to gei pe^ 


-3- 


one half 3 Recently, during a personal 
discussion, Th, von Kármán suggested to the author that the theory can 
be generalized to use the tangent at the state of pa corresponding 
to the undisturbed parallel flow. Thus the range of usefulness of 


Barad oan be greatly extended. This is carried out in the first 
arci, 


part of the present kantin 1. 
... ˙² this theory, based upon 
Demtpchenko and Busemann's work,is applied to the case of flow over 
elliptic cylinders and the resulte compared with those of 5,0. Hooker 
(Rot. 28) and C, Kaplan (Ref. DN Furthermore, results calculated by 
GlauertWa-Prandtle linear theory are also included for comparison. 
Reomtly, f. Bateman (Ref, A) demonstrated a remrkable 
reoiprooity ef two fields of flow of two fluids related by a certain 
point transformation, It will be shown in — this 
paper that the flow of incompressible fluid and the flow of compressible 


the 
fluid approximated by the use of tangent to adiabatic pressure-volume 


curve can be interpreted as such a point transformation, It is thus 


possible to obtain a solution for compressible flow whenever a solution 


of incompressible flow is know, — SEE, 
„„ i ^ 


incompressiblg flow can either be obtained by the well-known method of 
atisa) 


stios at 


me $ is the pressure, 1 is the specific volume and 


T ie the ratio of specific heats of a gas, the adiabatio relation 


pul constant is a curve inthe É-i plane as shom in 


Pig. la, Ta» conditions near the point ( h, i, which corresponds 
to tie stato of fiw undisturbed flow oan be approximated by the 
tangent to the curve at that point, The equation of the tangent at 
this point can be written as 


bo Clan)» cle 
ases gw! sem tha d.a ftoi — 


where U ia the slope of the tangent and (f de the density of the 
wer Tne slope ( must be equal to the slope of the curve at 
the point (J,, . therefore, 


0-02) - GHE)- LEIA 
where 4, is the See DEE the conditions 4, ,. 


Thus eq, (1) can be written as: 


f iid ste cS £ 
= nimaki T 
This is an 8 
is shom in Fig, E coe adiabatic relation, 


2 theorem for oompressible fluids is 


J 
nen n o 


$ 


where 44 is the velocity of the gas and the subscripts 2 and 3 
PA denote two different states of the fluid, By substituting eq. (2) 
into eq. (5), the following relation isdtained: 


Now it 4520, =u Big fac $ w 
the subseript O  denotée the state of|mest, eq. (4) gives: 


-s- 


If the square of sound velocity 5 


d ie defined ( 


usually done) as the derivative of Z with respect to ges 1 
eq. (2) gives: 


ag = As A efe 


Therefore, eq. (5) can be written as: 


9 
* /* = i 


ne 
It is interesting to m ha from eq. (8) the density 
Vedi catt 
(decreases as, !Teloctty increases, as EE Thus eq. (6) shows that 


Je att the local velocity of sound increas the velocity increases, This 


ad yan 


i» just opposite to the real gas, because in the oase of an adiabatic 
flow of a real gas it is well nom that the temperature of gus 
decreases as the ranri ihs — and thus the local 

sound volocätyJalso decreases, However, in the present approximate 
theory, the ratio Z or Mach's number still increases as the velocity 
increases as can be seen by eq. (7). But this ratio only reaches the 
value unity when g =O  ,orfro m. (8) wh e . Itis 
thas seen that the entire regime of flow is subsonic and thus the 
differential equation Vias —— This 
is the reason why the complex representation ot Ferootty potential and 


ju 


‘stream function is possible for all cases, as will be shown in Ae 
following paragraphs. However, one should realize that the portion of T 
tangent that could be used as an approximation to the true adiabatic 
relation is that portion which lies in the first quadrant, Thus the 

upper limitle velocity, for practical application of the theory dm OLD 


pro - By using eqs. (17) and (18), this upper limit 
is found to be 


Ehu ho LT „ 


Since the point (48) lies on the eerten: curve, the 
relation 74I can be used, and eq. (9) becomes: 


ei dutt 


it is not likely that the ratio ei will rise to values much 


RUN 


higher than ? 3 SE IQ MERCI LT 


where K. V are the componente of a abya 5 
respeotively, The equation of continuity, 


AU $k) ~ 


will be satisfied, if the stream function Y Ae introduced such 


that 
$ 2 ay 
f, “7 : -Í'— 


Now if the angle of inclination of the velocity 44 to 
the Y axis sf » eqs. (11) and (12) give: 
dé -w oop du + , ee, 
eee, e, 0m 
Solving for dy and 4 , 
e 2 E 
dp- SA. L. f. e, 
P au, a pdf 
So long as the correspondence between the physical e and hodograph 


plene v. one to ono, or mengen Ef fo. 
expres: and jas functions of ur ` dies? œa 0 


(aa) 


as functions of W, j . Thus, 


dh = du dw 4,48 (as) 
AH podot d, di > 
where primes indicate the derivative 8 — 


E) tastean 

indicated as subscripts, Now substitutiyeq. (15) into eq. (14), Pier de ‘4 

u al- et ae [P apnd (et 
HOLM EEA he (AG + pln 


Since the left-hand side of eqs. (19) are eyact differatisie, rs 
epaly the reciprocity relati. 


Ann) A IM EN) 
een) le) / 


Carrying out these differentiations and simplifying with 
the aid of eq. (7), eq. (17) gives: 


EHE we te ete Seb 
Käch PL = AU AH 


Solving for Pu 4’ ` 
2 


gis zu DN as) 
Eq. (19) can be further simplified by introducing a 


„ 28) 


new variable CÙ , such that 


40 Z EA ted 
Then eq. (19) becomes: 
d 


tho 
This can be easily recognized as E FienanibCauchy* differential 


equationsand tus +Å D must be an analytic function of 


(21) 


DER . However, for convenience of calculation, another new 
wet er independent vartaties Je W caß * 
y- W muß are introduced where J = ge 


Then eq. (21) can be written as: 


20 av) ia 


3) ~~ E 
— integrating eq. (20), 


W- Z en 


y au +4 


lade = 4M (24) 
> aw” ; 
Substituting into eq. (8), veten d, phe ban for de d satio = 
Rl dw ee 
S 4. W? (25) 


Eqs. (22), (23), (24) and (25) are the basic equations of 
the present theory. One recognizes eq. (22) as the RiemanntCauchy 
differential equationé, and thus the complex potential 7 = Ø 6, 
must be an analytt function of MW = Yu „or: 


puf- FIH-ıW)- Ffw) 

-if - Elan) ffn 
vor. it is necessary to find the wlues of E maf 
— eege Moss H . e.. 
to find the transformation from hodograph plane to physical plane. 

By using eqs. (24) and (25), eq. (14) can be written as: 


A: Se 2 2 22 


(27) 


El: zl ED fog 


2 
* w= UIV + These oquatione ean be combined into 


one equation by means of eq. (26). Thus, 


; ^ 
h = duh: . = es (28) 
Hence, if an analytic function Fr) is given for each value of 
W . the corresponding real velocity 44 can be calculated by 


eq. (24). 
which this velocity oocurs can be calculated by means of integrating 


Then the coordinate of the point in the physical plane at 


eq. (28). The pressure at thie point is then given by eq. (2). However, 


there is an objection to this procedure. That is, the invostigator 
a en 


nath te A 5 
ELE LE with the chosen function Z/W), we 
2 the desired shape ef selle boundary and flow pattern , 


———— In other words, this procedure still suffers 


the difficulty of boundary conditions as de is common to all hodograph 


methods, ay aa 


However, due to the particular simple relation of eq. (28), 


ome ascertainiapproximatel; e oi by 
starting with the el function, 


Flit)= dug - W Gl) 


(29) 


where M? is the transformed undisturbed velocity to be interpreted 


by eq. (25), and C, ie the complex coordinate fu . This 
fluid over 


function is so chosen as to give the flow of inoompres: 

the desired body shape in coordinates Fand 7 $ ä | 
r2 Verl 

Fhe real velocity 4 Ci plane of the incompressible fluid as 

transformed velocity ` Ä in the hodopraph plane fol/Gompressible 

fluid) 1. is knom that 


7 AGO 


| WM (50) 
S Thus w= N 2 


E t e eu 
wid, $ 


spear eqs. (30) and (31) se oq. (28), — 


Integrating, 


S u 
geb, & (se) 


Therefore, the complex coordinate in the physical plane 
ass fluid is equal to the corresponding complex coordinate \ 
in the physical plane ofVincompressible fluid plus a correction term. 

Since this correction term is usually small, the resulting shape of the 

body will be quite similer to the one inVfncompressible fluid. The 

factor|in the correction term depends upon the Mach's number of the 

undisturbed flow only. This can be pcg mus of eqs. (7), (8), 


and (23), because from those equations eme-hme the following relations Minus: 


RI 1 cn en 
= (My? _ 
OK ee fer (33) 


ka 72 D se ru mer of Mie tere flow. pA 222 Fa 


Mpi Glen d Nite. ici Kee planes I is 


first obtained from eq. (30) and then with eq. (23): 


PIE, 7 E AM (34) 
[7 ul ee, 


4 
Fow wer m Eiiphi Gylinders ae 
E theory dersiepal in-Pari=I will be 
applied to calculate the flow over an elliptic cylinder at zero angle 


of attack, The incompressible flow over an elliptic cylinder in the 


complex coordinate Č oan be obtained by applying Wwe Joukowsky's 


transformation to the flow vr a oiroular eylinder, in the 2 
Taa KC d 
coordinate S with ‘Located at the Therefore, 


— FM) or M, Ge can ve written as: | 


F-m Leis s) | 


4 = UT (38) 
Ee EAT 
„ on 


It is omvmimt to carry out the computation by using the 5” 
coordinates. Thus eq. (32) is rewritten in the following form: 


4 
" al s 
zz) [ a) p (88) 
c 
PU 
if Vd I ANE over 
the surface of the elliptic cylinder emi, are eegen, theu 


Z'- be? 
— 7 (39) 
as ot 
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where Ê is the argument as shown in Fig. Gand Ê is the radius of 
the oiroular section in the Z’- plane which determines the thickness 
ratio of the elliptic section in the  - plene, Substituting eqs. (36), 
(37) and (39) into eq. (38), and carrying out the integration, the 
following expressions for the Y and E coordinates corresponding to 

Z’ am er separating the real and imaginary parts, 


(^ ardent" Ay if bib pours I, eet 


Gs li- 7) Jed) An / HIT ES 744 { ar 


The horizontal and vertical semi-axis of the elliptic section can 


then be calculated by substituting f a ona J-F 


respectively into T (40). The thickness ratio d ie thus obtained 


ei 2 de MED . í 
1- ae, mag) 0m | 


For a given thickness ratio and Mach's number SEEN, adi 
n the value of NS, 1+ tiret computed by nonna of oq. (38), and ^ 
then eg. (41) is solved graphically for 4 » p 
After 4 is obtained, the coordinate 7 and fer 
each value of ÀJ can be computed by using eq. (40). It is fortunate 
that the values of X, 7 so obtained lies very close to the true 
elliptic seet4en 4... Hence, the velocity and the 
pressure distribution obtained by using eqs. (34) and (35) are considered 


/S6 

nn 
Im fl na mèh velina at Laver Inoki munter, UR 
uam Ate al dos a retell g TÉ eant 


If the velocity of flow over a body is gradually increased, 


the maximum local velocity in the field will also be increased. When 
2 
the maximum 


D 
local velooity reaches the local velocity of sound, Sen shock wives 


‘appear su the drag of the body suddenly ese, Danter: 


Pis velocity nu ‘considerable interest to practical engineers and 


is usually referred to as the critical velocity of the body. It is 
shown by Kaplan (| Ref. D and others that at this critical condition 


the mtio of maximum velocity Ay, in the field È RER 


velocity A , is related to WB Mach's number E Berl 
2 simaan in the following manner: 


4 
Cd. Er Hl (42) 


Wray, in the flow over an elliptic cylinder at zero angle of 
attack occurs at the top of the cylinder. Using eqs. (34) and (33) 


the value of “m is found to be: 
` 


EEN R t 
SET 
PATT RUM ay 


Equating eqs. (42) and (43), the equation for calculating 
the critical Yach's number — TZ er the undisturbed for 


sach value of $ ist 
( TET d 
SE 
KOSCH 
Lied 7 re" mi 


Knowing GL for each value of É , the 
corresponding value of J can be calculated by means of eqs, (34) 
and (41). Fig. e shows the result of this calculation with Kaplan's 


A Laer wale SS cation 


Mer 


value inoluded for "PT TT 
< EE a R 
A is oe a R Mis Ar NES 
offost of compressibility daa — 
consiatent with the rosita shom in rige. Odi ,, 
g 


Sehen 2 
The Mee of Ler ay Fan chons 


* bet aen Potion, E set F a ee ty 


kt e, 
A JV RSE y A hak ee 


eee bet Ad means of eqs, (11) and (12), it oan 
be show: that the following relations holdı 


| Ai: [pg sera (fe sully sd 
| AY reg [gere pad frin i 
4 m um that Lene, g Ki A ag p 


(7 9 wll eee, g bi. 


„„ f A em 


72 „ oe 

EN a 

is dd da, Ur 
KA AER 


j * A — dams , WIE , 
| RON r re 
L . 76 Ki D 
E ZA (M) 


Mae ai RL fn 4 A pis 


T2: d H-- fu 


Tram, The quantities Æ and „f have the dimension of a 
d 


c x 
velocity and se REN con be considered SAH velocity 


in the plane whose coordinates are denoted ty VY ana / . tis 


(49) 


relation between the ZL plane and the / amd La plane is shom 
ur 
in Fig. It is thus seen that if the undisturbed 
X Z plane is in the positive AT -direstion, the undisturbed t.. Mars 
in the Y -plane will be in the negative Me direction. Furthermore, 


f Soke thare exists a complebe reciprocity habısan the Ki plane and | 
L Í the Y Y piani thom by H. Batemın (wr. 24). 


ee at w t ioka aai X Fans. 
P 
GE E 


apply to fluiés of erbitrary properties, However, since only the flow 
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6 wll Ann 

of incompressible fluid „ At would be 

Saane 9 to G 771 the YY plene it 
A 

the Sex in the 4 y plane is incompressible. If the fluid in the 


E ^f plane is incompressible, then = =/ ‚end the Bernovllitg C 
ete: Aer. gies: 


+ = = 


(82) 

Wem Ts ine Y e et 5 — . 
7 3 then Mernoul1I em gives: 

e SE (ss) 


Ma EF 


—— In view of ede. (49), (50), and (52), eq. (53) oan be 
written in the following form: 


+ 
re) Al SS = constant (54) 
2 5 By differentiating eq. (54) with respect to 5 D 


multiplying the — expression by = and then integrating 
with respect to SC .— the following relation conneoting 


the pressure P and Rensity C~ for the fluid in the VY planes A. 


P= ml. ife. (55) 
whee C the wë SÉ 


~e Comparing eq. (55) with the approximate adiabatic relation 
eq. (2), also noting eq. Ya). it is evident that eqs. (55) and (2) 


are identical if: 


bi 
© 


LEA L, GST 
se E. Rit het 


In eq. (56) Ä As the velocity of sound in the H-V plane, and the 
"E 1 refers to the conditions in the undisturbed pad 

ci 
Henco, 4 is wa Mach's number of the undisturbed 


ü 
— — Ano 


j e 
N Jar mi + Q, thy thiu be eee, Zen gp C) 4 (4) 
Bb ts SW 


D 


LE ur RI > mA 
th ar 
p 4 (St) em eE, aa: 


(74 
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(T 7- KÉ Ee cof atat a, 


8: bit f T F , 
E, 4 a pu Té fg P 757 at ët 


uug ukla s 
ly ~ WK 
TVA A 
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OF COMPRESSIBLE FLUIDS 
by 


Heue-shen Tsien 
California Institute of Technology 


Introduction 


Assuming that the pressure is a single-valuod function of 
density only, the equations of two-dimensional irrotational motion 

of compressible fluids can be reduced to a single non-linear equation 
of the velooity potential. In the supersonic case, that is, in 

the oase when the flow velocity is everywhere greater than that of 
lool sound velocity, the problem is solved by Meyer ^ Prandtl and 
Dusemann using the method of characteristics, The 


ential diffioulty 
of this problem lies in the subsonic oase, that is, in the case when 
the flow velocity is everywhere smaller than but near to the local 
sound velocity, because then the method of characteristics cannot 
be used, Glauert + Prandtl (Ref. 1) treated the case when the 
disturbance to the parallel rectilinear flow due to presence of a 
solid body is small. They were then able to jinearize the differential 
equation for the velocity potential and obtained an equation very 
similar to those for the incompressible fluids. But there are usually 
stagnation pointe either in the surface of the body or in the field 
of flow, where the disturbence is no longer mall. "enge, it is 
doubtful whether the linear theory can be applied to the flow near a 
stagnation point, Om the same ground, the theory breaks down in 
case of bodies whose dimension across the stream is not small compared 
with the dimension parallel to the stream, 

To treat casos in which the body is blunt nosed, Janzen 


and Rayleigh developed the method of successive approximations. This 


method was explained physically and put into a more convenient form 


by L. Poggi (Ref. 2) and P. A. Walther (Ref. 3). Recently C. Kaplan 


(Ref. 4) treated the case of flow over Joukowsiy airfoils and elliptic 
cylinders, using Poggi's method. However, the method is rather tedious 
and the convergent very slow if the looal velocity of sound is approached. 

Molenbroek (Ref. 5) and Tsehapligin (Ref. 6) suggested 
the use of the magnitude of velooity Ai and inclination £ of 
velooity a chosen axis as independent variables, and wore thus able 
to reduce the equation of velocity potential to a linear equation. 

This equation was solved by Tschapligin (Ref. 6) and recently put 
into a more oonvenient form by F. Clauser and M, Clauser (Ref. 7). 
The solution is essentially & series each term of which is a produot 
of a hypergeometrio function of déi and a trigonometrio function 
of f + The min difficulty in practical application of this 
solution is to obtain a proper set of boundary conditions in the 
plane of independent variables W, A and put the solution in a 
closed form. 

Tachapligin (Ref. 6) shows that a great simplification of 
the equation in the hodograph plane results if the ratio of the specific 
heats of the gas is equal to -1 . Since all real gases have their 
ratio of specific heats between l and 2, the value -l seems without 
praotical significance. It was Demtschenko (Ref. 8) and Busemann 
(Ref, o) who made olear the meaning of this specific value of -l, 

They found that this really means to take the tangent of pressure- 
volume curve as an approximate to the ourve itself. However, they 
limit themselves to use the tangent at the state of rest of tho gas, 


Thus their theory oan only be applied to flow with velocitit 


up to 


one half that of sound velocity, Recently, during a personal 
discussion, Th, von Karman suggested to the author that the theory oan 
de generalized to use the tangent at the state of gas corresponding 


to the undisturbed parallel flow. 


Thus the range of usefulness of 
the theory oan be greatly extended, This is carried out in the first 
part of the present paper. 
In the second part of this paper this theory based upon 
Dembschenko and Busemann's work is applied to the oase of flow over 
elliptic cylinders and the results compared with those of 5.0. Hooker 
(Ref.10) and C. Kaplan (Ref. 4). Furthermore, results calculated by 
lauert e-Pren at! % linear theory are also inoluded for comparison, 
Recently, H. Bateman (Ref. 11) demonstrated a remrkable - 
reoiprooity of two fields of flow of two fluids related by a certain 
d point transformation, It will be shown in the third part of this 
paper that the flow of incompressible fluid end the flow of compressible. 
fluid approximated by the use of tangent to adiabatio pressure-volume 
curve oan be interpreted ae such a point transformation, It is thus 
possible to obtain a solution for compressible flow whenever a solution 
of incompressible flow is known. The difficulty, however, lies in 
the proper choice of this incompressible flow so that the solution of 
i compressible flow will have a solid bowidary desired. Fortimately, this 
i transformed boundary is independent of the Mach's number of undisturbed 
` flow, Therefore, by a modification of the procedure, two relatione 
between the velocity and the pressure distributions over a body at low p 
1 velocities and those at high velocities ver the same body are obtained, 
This enables one to predict the oharacteristios of the flow over a 


ed characteristics of 


body at any Mach's number whenever the low 


the flow over the same body are known. The characteristics of the 


incompressible flow oan either be obtained by the well-imown method of 
conformal mapping or by experiments. 
merodynamio enginecrs usually have the low speed charasteristios at 
hand and that high speed data have to be obtained by use of a costly 


Due to the fact that practical 


high speed wind tunnel the above mentioned relations are believed to 
bo of considerable use to them, 
In the fourth part of the paper, the theory developed in 
Part III is applied to correlate airfoil data obtained by J. Stack (Ref, 12) 
in the 5.A,C.A, 24" high speed wind tunnel. The agreement with theory 

is found to be satisfactory, Then this theory is applied to predict 


the comprossibility effect on the lift and moment of N. A. b. A. 4412 
airfoil using experimentally determined pressure distribution over 
the same airfoil at low speed. The result is again compared with the 
more simple Glauert-Prandtl theory. 


1, L 


if f is the pressure, / is the specific volume and 
is the ratio of speoifio heats of a gas, the adiabatic relation 
yle constant is a ourve in the fr plane as shown in 
Fir, la. Now conditions near the point (h, v) which corresponds 
to the state of fluid at undisturbed flow oan be approximated by the 
tangent to the curve at that point. The equation of the tangent at 
this point can be written as 


p-p- Giele cls 


where (7 is the slope of the tangent and © ie the density of the 
gas. Now the slope ( must be equal to the slope of the curve at 
the point f 71. — , therefore, 


c- (E) tw 
were d isis red ef E X 9s 


Thus eq. (1) oan be written as 
CEP 
DN LEE x 


This le an approximate pressure-density to adiabatic relation, and 
is shown in Fig. lb with true adiatatio relation, 
The Bernoullis! theorem for compressible fluids is 


2 
J l- [ 4 © 
iss, <5 


where Ad is the velocity of the gas and the subscripte 2 and 8 


denote two different states of the fluid. By substituting eq. (2) 
into eq. (3), the following relation ischtained: 

a LL eb 
A RU. cu rey os £ 

How if wé, r. yo f, mä = 
the subscript © denotes the state of rest, eq. (4) gives: 
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If the square of sound velocity Æ is defined (as 


usually done) as tho derivative of D with respect to © , 


eq. (2) gives: 


2 4s = e$ constant 


Therefore, eq. (5) can be written as: 
ues EN 
GI -/- $ in 


Similarly, (8) 


7. 
[S Ar 
It is interesting to notice that from eq. (8) the density 


deoreases as velocity increases, as expected. Thus eq. (6) shows that 


the local velocity of sound inorence 


s the velooity increases. Thie 


1s just opposite to che real gas, because in the case of an adiabatic 
flow of a real pas it is well known that the temperature of pas 
decreases as the velocity of gas is increased, and thus tho local 


sound velocity also decreases, Howover, in the present approximate E 
Ja 7 theory, the ratio a Or Mach's number still inorcases as tho velocity 
2 increases, as can be seen by eq. (7). Sut this ratio only reaches the 
dab à, tn value unity when ¢ =  , or from eq. (8) when “7 e, It ia 
. thus soen that tho entire regime of flow is subsonic and thus the 
2 E v differenti&l equation of velocity rotell. This ' 
P is the reason why the complex representation of velocity potential and i 
Za d stroan function is possible fcr all cases, 4 will be shown in 
Ka 7 hy TSL Lein paragraphs. Towover, one should realize that the portion of | 
f fon ` tangent that ceuld be used as an approx ration to the true adiabatio 


relation is that portion which lies in the first quadrant. Thus the 


upper limit of velooity for practical application of the theory 1a 


when 7 =0 . By using eqs. (17) and (18), this upper linit 
is found to be 


(he Hkh „ 


Sinoe the point ., J, ues on the true adiabatic curve, the 
relation 4, = rt can be used, and eq. (9) becomes: 
7 


— 
C2» A Gu 22 E (10) 
This relation is plotted in Fig. 2, Since for most practical oases 
4t le not Likely that the ratio GI will riso to values moh 
higher than 2 D will remain positive, and thie theory will give 
an approximate solution. 
If the flow is irrotational, there exists a velocity 
potential $ such that 
&, E ar Q1) 
where L, V are the components of “ in Y Se 
respectively. The equation of continuity, 


Di YT D 


will be satisfied, if the stream function 4 is introduced such 
that 


2 
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Now ir the angle of inclination of the velocity “ to 
the X axis is d » eqs, (11) and (12) give: 
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Solving for fy, and i 
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(14) 


So long as the correspondence between the physical plene and hodograph 


plane is one to one, or mathematically 32H Ho > one can 


as funotions of 4 ‘A and so fu 


express y o 


as funotions of 4. . Thus, 


d4 = 2 du + 44 (15) 
df- Yor toa + Yy dt 
where primes indicate the derivative with respect to variables 
indioated as subscripts. Now substi tutined. (15) into eq. (14), one has: 


bs (Ld MEA ae tl Ye Séis 
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Since the left-hand side of eqs.(16) are exact differmtials, one can 


apply the reciprocity relation, and obtain: 
ilz ei 20e, a, run, 

cop * „dan 
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Carrying out these differentiations and simplifying with 
the aid of eq. (7), eq. (17) gives: 


. em $n HH E. EN 
Sb, nde 7 (18) 
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Solving for d. anè M 
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Eq. (19) can be further simplified by introducing a 
i new variable % , much that 


L du=£ 4 e í 
j Then oq. (19) becomes+ GER DN 
e > 

te ti 


This oan be easily recognised as the Riemann-Cauchy's differential 
equations and thus jf mast be an unalytio function of 
(D-i . However, for conveniense of calculation, another new 
set of Andependent variables //= JV anf " 
V - W oA are introduced where Wane” 


Then eq. (21) oan be written as: 


4 - xp 
| db. 


and also by integrating eq. (20), 


(22) 


Alo w = 
W = Je p. (23) 
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EV 


H "3 
U R w: 
Substituting into eq. (8), we have: 
& u 4W* 
P x hep = (25) 


Eqs. (22), (23), (24) and (25) are the tasic equations of 


4 the prosent theory. One recognizes eqs. (22) as the Riomann-Cauohy 
differential equations, and thus the oomplex potential = A. 74 
mé de aids Nesse o AP „or 


grip- rp(t)- Hm) 

4-il- Ef Mn) = Fin) 

Now, it de necessary to find the wlues of V and 
corresponds to a given set of values of 7/ and , , ies, 


to find the transformation from hodograph plane to physical plane, 
By using eqs. (24) and (26), eq. (14) oan be written as: 


er SE E | 
u^ 2 ent J. 2 2 


where we My . These equations oan be combined into 
one equation by means of eq. (26), Thus, 


den An = = H 


Zar (28) 


Henoe, if an analytic fmetion An is given for each value of 
wW , the corresponding real velocity +t” oan be oalculated by 


(26) 


eq. (24). Then the coordinate of the point in the physical plane at 
which this velooity occurs can be calculated by means of integrating 
eq. (28). The pressure at this point is then given by eq. (2). However, 
there is an objeotion to this prooedure, That is, the investigator 
has no way of telling whether with the chosen funotion Hi) . he 
is going to obtain the desired shape of solid boundary and flow pattorn 
which he is interested in, In other words, this procedure still suffers 
the diffioulty of boundary conditions as it is common to all hodorraph 
methods. 

However, due to the particular simple relation of eq. (28), 
one oan ascertain approximately the resulting shape of the body by 
starting with the particular function, 


Hi) = $n ZAC 66) (29) 


where M; is the transformed undisturbed velocity to be interpreted 
ty eq. (25), and Z ie the complex coordinate Fry . This 
function is so chosen as to give the flow of incompressible fluid over 
the desired body shape in coordinates F and + Now interpret 
the real velocity in [-/ plane of the inecmprossitle fluid as the 
transformed velocity ` 20 in the hodograph plane for compressible 

fluid, Then it is known that 


77 4 6C 
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Substituting eqs. (30) and (31) into eq. (28), one has: 


ae) 72 
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Integrating, 
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2 6-2 0D) LEE (32) 
Therefore, the complex coordinate in the physical plane 
of compressible fluid 16 equal to the corresponding complex coordinate 
An the physical plane of incompressible fluid plus a correction term. 
Since this oorreotion term is usually small, the resulting shape of the 
body will be quite similar to the one in incompressible fluid. The 
factor in the correotion term depends upon the Mach's number of the 


undisturbed flow only, This oan be seen by means of eqs, (7), (8), 
and (23), because from those equations one has the following relation: 


Ly. 2 
KA? fr 277 (33) 
whero 15 ia the Mach's number of the undisturbed flow. 
To oaloulate the velocity in physical plane, Zi 1s 


first obtained from eq. (30) and then with eq. (23): 


ae Jidd —E, = (34) 
EI Oe: " 
ET 
Then the pressure can be onloulated by using eq. (2). 


With some calculation the following expression for pressure at any 
point is obtained: 


IE og [TT e 


In this part, the theory developed in Part I will be 
applied to calculate the flow over an elliptic cylinder at rero angle 
of attack. The incompressible flow over an elliptic cylinder in the 
comolex coordinate S, can be obtained by applying the Joukowaky's 
transformation to the flow over a olroular cylinder in the complex 


coordinate Z with ite center located at the origin. Therefore, 


tno tuotin FM) or AGE) oan ve miren ae: 


(37) 


It is omvenient to carry out the computation by using the er 
coordinates, Thus eq. (32) is rewritten in the following form: 


1b) 2% (der de 
2= (2 EIN E EA (38) 
d 
If one limits one’s self to oaloulate the conditions over 
the surface of the elliptic cylinder only, 


Z- bet 
g= be? 


where Â 1a the argument as show in Fig, Sand 4 is the radius of 
the circular section in the g plane which determines the thicknoss 
ratio of the elliptic sootion in the Z - plane, Substituting eqs. (36), 


(37) and (39) into eq. (38), and carrying out the integration, the j 
following expressions for the I and Z coordinates corresponding to 


4 ST are obtained by separating the real and imaginary parts, 
d SC ZA, Ae Dä 
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The horisontel and vertical semi-exis of the near elliptic section oan 
then be caloulated by substituting 0 ma I-E 
respectively into eq. (40). The thickness ratio 4 is thus obtained 


1. HE eee ed I+ Zéi 3 s e e " (a1) 
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For a given thickness ratio and Maoh's number of undisturbed stream, 


the value of + Gy? is first computed by means of eq. (35), and 
then eq. (41) ia solved graphically for 4 


After Á ie obtained, the coordinate I and E: for 
each value of / oan be computed by using eq. (40). It is fortuna te 
that the values of 4, so obtained lies very close to the true 
elliptic section as shown in Fig, 4 Henco, the velocity and the 
pressure distribution obtained by using eqs, (34) and (35) are considered 


as those over the true elliptic sections. 
ration, dor 
n wa (ei „ are carried out and the results shom in Fig. 6, 


) Calculations for two thiaknı 


Fig. 6, Fig. 7 and Fig. 8, together with those of Kaplan (Ref. 4) and 
Hooker (Ref. 10) included, In order to compare the resulta with that | 
obtained from less exact theory of Glauert and Prandtl, these compu- 
tations are made (Ref, ), and the results also included in Fig. 5 
to Pig, 8, It is sem from those figures that the effeot of compressibility 
of fluids shown by the present theory is more pronounced than those | 
1 shown by Kaplan's and Hooker's calculation. In view of the fact that 
j Kaplan and Hooker carried their oaloulationmly to first approximation, 
and higher approximations will inorease the walue somewhat, it is 
probable that the walues given by the present theory are nearer to the 
true value, As to the results given by Prandtl-Glauert theory, it is 
seen that the theory breaks dom near the stagnation point and the 
results are poorer for thiok sections than those for thin sections, 
However, in oase of thin sections this very approximate theory does 
give satisfactory information for a point not too olose to the stagnation 
point. 
If the velocity of flow over a body is gradually increased, 
the maximum local velocity in the field will also be increased. chen 
a certain velocity of the undisturbed stream is reached the maximum 
local velocity reaches the local velocity of sound. Then shook waves 
usually appear and the drag of the body suddenly increases, Therefore, 
this velocity is of considerable interest to practical engineers and 
is usually referred to as the critical velooity of the body, It is 
shown by Kaplan ( Ref. 4) and others that at this oritical condition 
the ratio of maximum velocity far in the field to the undisturbod 


velooity ^ » is related to the Mach's number 4 


stream in the following manner: 


Ga lire Hef" (42) 


nay in the flow over an elliptic cylinder at zero angle of 
attack ours at the top of the cylinder. Using ede. (34) and (35) 
the value of 2 45 roma to bei 


(Gy 
T. LEJI = 
E77 22777) 7- ray 


Equating eqs. (42) and (43), the equation for mloulating 
the oritioal Mass number (. ge of the undieturbed stroan for 


mh wine et 4 im 
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Al rtt emp d. 


Ynowing bui for each value of 4 . the 
corresponding value of f oan be calculated by means of eqs. (34) 
and (41). Fig. 9 shows the result of this oaloulation with Kaplan's 


zë Ze 47 . fer (49) 


Therefore, tho quantities E and L have the dimension of a 
velooity and so R and S can be considered as a new velocity 


in the plane whose coordinates are denoted ty YY ana / . This 

relation between the 24 plane and the amt plane is shom 

in Fig, 10, It is thus som that if the undisturbed stream in the 
XY plane is in the positive E- direction, the undisturbed stroam 

in the Y . plane will te in the negative Kaeeetie, Furthermore, 


if KAS bE 


— 
% 


At sur? 


one has from eq. (48), 


So far the relations obtained are general, i.c., they 
apply to fluids of arbitrary properties, However, since only the flow 


of incompressible fluid has been studied exhaustively, it would be 
interesting to find out what will be the flow in the I H plane if 
the flow in the Ri plane is incompressible, If the fluid in the 

X $4 plane is incompressible, then $-/ s and the Bernoulli's 


H- $2 x d (52) 


Words Bee YY tmi, 16 77 ambe Mia rear and e denote 


2, 1 then Somoulli's theorem rives: 


JE +H QI e (ss) 


In view of eqs. (49), (50), and (52), eq. (53) oan be 
written in the following form: 


4/4) dP T d 
SR TES A PE ES 
By differentiating eq. (54) with respect to Z D 
multiplying the Dees expression by T and then integrating 


— — 


with respeot to = a une tas the following relation comesting 

the pressure P and density 9. for the fluid in the Y I, plano: 
— 

P= omens -4 220 = (55) 


Comparing eq. (55) with the approximate adiabatic relation 
eq. (2), also noting eq. (16), it is evident that eqs, (55) and (2) 


are identioal if 


LEA ALY] 


> 
and Constant = Rita 


In eq. (56) 4 is the velocity of sound in the / J, plane, and the 
subscript 1 refers to the conditions in the undisturbed parallel 
stream, Hence, 


is the Mach's number of the undisturbed 


i parallel stre. ` ` 

N To find the coordinates Y and Y in terms of 

T LA . eq. (48) has to be integrated, It is convenient in this oase 

to use the complex potential of Feste flow in the I pimo. - 
Thus, if 


dup- ba- ar Gl2) 


(57) 


j Then it oan be easily shown with the aid of eq. (52) that 


Fa db m 


It le interesting to notios that eq. (58) is identical to the 
correction terms in eq. (82) except a factor, Tiowevar, there is a 
fundamental difference between the method explained in Part I and the 

method explained in this Part, The previous method will give different = 
solid boundaries for different Mach's number of the undisturbed flow 

even if the samo function (G 18 used, Eq. (58) shows that if the 


seme function (7 1s used, the resulting solid boundaries are all 
geometrically similar. 


By using eqs. (52) and (49), the componente of velocity 
in the JN plane can be expressed as: 


In eqs. (59) and (60), it is assumed that the undisturbed flow in 


the V plane is in the y-direction and so the undisturbed flow 
in the’ J, plane is in the P ätrsotion. The relation between 
PETS DE then be obtained from eqs. (56) and (57), that io: 


the EE SE - Hf) 


(01) 
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Thus, eq. (60) oan be rewritten as: | 


2 i 
Ep (62) 


For à given oddini of fler aa te E planos eqs: (88) end (61) 
will determine the characteristics of the flow in the Z plane. 
Using eqs. (55), 55) and (50), the pressure in the ff Y plano oan te 


Thus by means of eq, (61), 


J- (HOT Zr 2 5s 


(63) 


Thus as the mothod is developed one starts with a certain funotion 
Gl2) roprosenting the Inoamprosstblo flow in the Af piano, we 
first calculate the coordinates . J, and thus the solid boundary 
in the E H plane. Then with the desired value of for flow 
in tho YY plane, tho velocity and pressure distribution can be 
saloulated by means of ede. (61), (62) and (63). Thus the flow problem 
of approximately adiabatic compressible fluid is solved by means of 
an auxiliary function G(2) for flow of incompressible fluids. 
However, this procedure suffers the defect of diffioulty in obtaining 
the desired solid boundary compared with the first method developed 
in Part I, because now Z is not related to 7 ty a small correction 
as is in oq. (32), but by a transformation given by eq. (65). 
However, a mich more useful form of eqs. (62) and (63) 

can be obtained Ar one remembers that the solid boundary of the flow 
im the F Y plano is independent of Iach's number, but only on 

Gle) . er in other words, on e . Now substitute oq. (61) 


into eq. (62) and let HI +e 
8 TT 


= *- od (64) 
wore [E] denotes ‚che valosity for flow in tne E plane 
when 7 t. Then substitute eq. (64) into eq. (62). The 
following relation between the velocity of incompressible flow and 
the velocity of compressible flow over the same solid boundary is 

obtained: 


(65) 


E 
where £4 is a function of CG) given by eq. (61). Similarly 
eq. (63) oan be written as: 


—7 
T= (+4) SA (ac) 
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1 Jf tor EE is denoted jy I, . then 1t ta well 
imom that 
r= -E: 
" Git 


Substituting eq. (67) into eq. (66), we have: 


-LAr 
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Eqs. (65) and (68) together with eq. (61) thus enable 
one to calculate the velooity and pressure distribution over a body 
in compressible fluid once the velocity and pressure distribution over 
‘the same body in incompressible fluid is mown. The latter oan be 
obtained by either well-known analytical means or by experiments. 

1f eq. (61) is substituted into eq. (68) and only first order 
terns of LÉI is retained, then 

d 


(69) 
Thus, 47 (t and 4 are the lift coefficients of an airfoil 


at (#) and at 4) = and the same angle of 


attack respeotivoly, one has from eq. (69): 


i = DE (70) 
Wr 


which checke with the Prandtl-Glauert theory (Ref. 1). 


T. 
By inverting eq. (68), one has: 
(m) 


1f the theory developed in Part III is true, then the value of I 


2 
‘ferent values of Mach's number t } of 
calculated fron // at aif: values o EP) 
undisturbed flow should fall on the same ourve, This is done for the 


case of upper surface of NACA 4412 airfoil tested by J, Stack (Ref. 2) 
at % angle of attack, The Jf ^re replotted in Fig. 10, while 
the M, curves onloulated from J curves are show in Fig. 11. 

The derived values of I, lie fairly olose to me ourve, Since the 


data for this calculation is taken from a small soale figure, the 
deviations must be largely due to inaccuracies involved. Tip. 11 can 
thus be taken as an experimental check of the theory. 
One of the immediate possible applications of the theory 
developod in Part III is to determine the oritioal speed of a body. 
If tfh mu nd Mma Aro the maximum speed und) 

C minimum pressure over a body, Then by eqs. (65), 
(87) and (42), the equation for oritioal Mach's number is obtained as: 
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or 
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In both eqs. (72) and (73), £ is a funotion of A given 
h by eq. (61). Eq. (78) is plotted in Pig. 32, together with a ourve 

given by B. k. Jacobs (Ref. 13) based upon/Glauert-Prandtl theory. It 
le thus seen that Jacobs's ourve rives a critical Mach's number, which 
is a little too high, This difference is, of course, to be expocted, 


due to the higher order approximation of the present theory. 
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SIMILARITY LAWS OP HYPHRSONIC FLOWS 
Í by Hsue-shen Teien 


Introduotion 


Hypersonic flows are flow fields where the fluid velocity is 
much larger than the velooity of propagation of small disturbances, 
the velooity of sound. Th. von Kärmen (Ref. 1) has pointed out that | 
A in many ways the dynamiı | 


of hypersonic flows is similar to Newton's 
corpusailar theory of aerodynamics, The pressure acting on an ine 


clined surface is thus greater than the free stream pressure by a 


quantity which is approximately proportional to the square of the 

N angle of inclination instead of the usual linear law for conventional 
supersonic flows, E, Sanger (Fer. 2) has, in fect, used this concept 
to design the optimum wing and body shapes for hypersonio flight at 
extreme speeds. 

Recently, von Karman (Ref. 3) has obtained the similarity laws 
for transonic flows where the fluid velocity is very near to the 
velocity of sound, He deduced these laws by using an affine trans- 
formation of the fluid field so that the differential equations of 
the flows are reduced to a single non-dimensional equation. In 

I this paper, the same method is used to derive the similarity laws 
] for hypersonic flows. These laws will be, perhaps, useful in oorre- 
lating the experimental data to be obtained in the near future by 


hypersonio wind tunnels now under construction. 


1 Differential Equation for Hypersonic Flows 


If u, v aro the components of velocity in the x, y directions 


end a is the local velocity of sound, the differentiel equations for 


irrotational two-dimensional motion are 
BIBER) +l ae 


Now if a slender body is present in an otherwise uniform strom of 
velocity Y in the x-direotion, equations (2) is satisfied by intro- 
duoiag the disturbance velocity potential y defined as 

a= Vt # 

v= X 
If a, is the velocity of sound for the gas at rest and a is the 


(3) 


welooity of sound corresponding to the free strom velocity V, 


then there ore the following relations: 
a " B 
ve BE LEHE o 
a 2 747 A 
a % A 
where y ie the ratio of the specific heate, 
For hypersonic flows over a slender body, both a? and # EN 
aro small in comparison with V, By substituting equations (3) and 


(4) into equation (1) and retaining terms up to second order, one has 
[ir min zë . PRSE: - uL 
+ [rye GE et «6 


Here M? is the Mach number of the free stream or 
17. A (6) 


a 


(8) 


— — adl cd 


Similarity laws in Two-D; l Flow 
lf 2b is the length or chord of the body and 4 the thickness 


of the body, the non-dimensional coordinates E] and 7 oan be de- 
d KMD 
n ral 


where n is the exponent yet to be determined, von Kármán (Hef, 1) 


fined as 


(7) 


has shown that for hypersonic flow over a slender body the variation 
of fluid velocity due to the presence of the body is limited within 
x narrow region olose to the body, the hypersonic boundary layer. 
Therefore, in order to investigate this velocity variation, one must 
expand the coordinate normal to the surface of the body. This is 
similar to the case of ordinary viscous boundary layer, where Prandtl's 
simplified boundary layer equation is obtained from the exact Havior- 
Stokes equations by a coordinate expansion normal to the surface 
of the body. From this reasoning then, n must be positive so that J 
is much greater than Go. This surmise is substentiated by the 
lator calculations to be shown presently. 

The appropriate non-dimensional form for the velocity potential 

q ie 
g- abi He v) [ 


By substituting equations (7) and (8) into equation (5), one has 


E D-rg- gp . e A 555 E 00 


pb Bt RET) . 


The boundary oonditions at infinity require that the flow 


velocity be W Thus 
# 5 D zo A Dei 


If the slender body is a symmetrical one, then the condition at the 
surface of the body oan be written as 


ap... = Mt) 165) m 
where MIT for -1« f < 1 isa given function describing the thick- 
ness distribution along the length of the body. Equation (11) oan 
be converted into the following form by means of equations (7) and 


(8). a porn 
GIL: wu") 2 (az) 


Since the body le thin, d) is very mall. Therefore the first 
group of terms in equation (9) is negligible in comparison with the 
rest, Then both the differential equation and the boundary conditions 
oan be made to contain only a single parameter if one sets 

nel Dal 
That is, ir 


WE Kk 
then equation (9) becomes 


% ette 


end the boundary conditions become 
e at 
ZZ: 
85 = KI om -14541 
ye 


The meaning of this similarity law is the following: If a series 
of bodies having the seme thickness distribution but different thiok- | 
ness ratios HA) are put into flows of different Mach numbers M 
NO such that tho products of MO and (4/j) remain constant and equal e 
to K, then the flow patterns are similar in the sense that they are 

governed by the same function //f,7) , determined by equations (14) 

and (15). 

If po is the stagnation pressure, p? the free stream pressure, 


d and p the looal pressure, then 


p-p Lert] ^ 
IE 


In notations introduced previously and by retaining terms of proper 
magnitudes, one oan write the expression for the local pressure as 


L ) 
Perl d Agee] De 
The drag D of the body oan then be oaloulated, It is given by the 


following Sura 


D- Bt IK]; = uf Wfl [A n-o vel M 


lf one wishes to compute tbe drag coefficient Cj, then one de 


7 
E. J Géi D -roi ^ Kae] imd T an 


PIE A 
dE a) wi 


For a given thickness distribution, the quantity within the braokets 


Therefore, one 


is only a function of K, the similarity parameter, 


oan write 


wi 
Similarly, one obtains for the lift coefficient Ki of the lift L 


5a 25 A(K)* ger. Med Q8) 


the following lawı 
l C. E E AC. pans) Qs) 
These similarity laws show that for bodies of the same thickness 
distribution at angles of attack proportional to the thickness ratio 
(4/b ), the quantities (ce) and SÉ) are functions of the 
single paramotor K = w0 L , 

Equations (18) and (19) agree with the results of the moro 
limited linearized theory of Ackeret (et, 4). Aooording to thia 
theory, for similar bodies in tho sense stated above the drag 
coefficient and the lift coefficient are given by 


WI 
For hypersonic flows of very large values of MP, these expressions 
reduce to D 
ar höhe (20) 


s fnr (21) 


Equations (20) and (21) agree with equations (18) and (19). Equa- 
tions (18) and (19) are, however, more general and complete. 


For axially symmetrical flows, the ordinate y is the radial 
distance from the axis to the point concerned, Then a similar 


analysis leads to the following differential equation and boundary 
conditions! 
H 


fE: Bea erf SBOE E 


gue d 
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were MË 7 is the distribution funotion for aross-seotional aroas 
long the length of the body and EN as for two-dimensional 
flow. 


(24) 


The drag coefficient Cp referred to the maximum cross section 
of the body is then governed by the following similarity lew: 


(28) 
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THE "LIMITING LINE" IN MIXED SUBSONIC AND SUPERSONIC rat 
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It is well-known that the vorticity for any fluid element is constant 


if the fluid is non-viscous and the change of states of the fluid is isentropic, 
When a solid body is placed in a uniform stream, the flow — of the 
body is irrotational, Then if the flow is further assumed to be isentropic, 
the vorticity will Pe sere en the while tha of flow. In other 
words, the flow is irrotational, For auch flow over a solid body ef finite 
a, it is shown by T. Theodorsen (Ref, 1) that the solid body exper- 
dences no resistance. If the fluid has a small viscosity, its effect will 

be Limited in the boundary layer over the solid body and the body will have 

a drag due to the skin friction. This type of essentially isentropic 
irrotaticnal flow is generally observed for a stream-lined body placed in 


a uniform stream, if the velocity of the wtream'is kept below the so-called 


" or rather at a certain value of the ratio of Z6 ken? 


Ge of sound, shock m appear, This phenomenw is called the 


“compressibility bufble' wr ro epe the change of states of the 

fluid is no more isentropio, although still adiabetio, Thés results in an 
increase in entropy of the fluid and generally introduces vorticity in an 
éripinalky irrotational flow. The inorease in entropy of the fluid is, of 
course, the consequence of changing part of the mechanical energy into heat 
energy. In other worüs, the part of fluid fffected by the shock mvo has 

much less mechanical energy. Therefore, with the appearance of shock waves, 
the wake of the stream line body is very much ea, and the drag increases 
drastically. Furthermore, the accompanied ehange in the pressure distribution 


E m 
over the body changes ent acting im it and in case of an airfoil 


decreases the lift force, 


H 
l.S L and the width of the square test section be b. Thes the Reynolds num- 1 
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he taf wb ds dm 


of extremely large viscofs effect. In other vords, the boundary layer 


W111 be so thick as to occupy the main portion ef the nos:le passage. To es 


demonstrate tnis effect, let we-esasiées het the length ef the test section 


ber based uon the conditions in the test sestion is Re = FE where U is 
the velocity in the test section. If, as a rough estimate, ve take the 
thickness of the boundary layer to be sero at the beginning of the test 
section and equal to a value | calewlated by the vall-known Blasius formule 
for » flat plate at the end of the test section, then 


tein o 


then 
2 * 
Tage „ KUE [C] 
Ou the other hand, tbe ratio of the meen free path Í amd the boundary 


X 
layer thickness ] 1s komm (Ref. 2) to be equal to ie 
Md 


Ze KE [D] 


where y is the reiio of sjesifio heats and cam be taken as l4, M io the 
Mach mmber in the test section. Ry combining (2) and (3), we have 
(lu olan u (e 

This relation is shown ia Fig. l. Thus for s Mach mmber M equal to 2, 
and L/b = 2, the boundary layer will completely fill up the test section, 
1f the mean free path is equal to 3.66 of the Youndary layer tbiokness or 
2.08 ef the tunnal width, This means thet the exiremely strong viscous 
effect at low derivities makes the ordinary sonsest of desigaing » vind 
tunnel totaliy inapplisehler 

Tos seme, Casp qon be also demonstrated by elu the ratio of 
frighiowa? loss on the wells ef the test section and the besk loss in the 
diffuser after the test sestien. rr 


tube of approximately the same cross-sesticmal ares as the test section, 
‘then tbe pressure lose dus te friction Af, is 


aur quern 


mim Ha eth h e techn k b 4 he uto f. sai, abe 
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Now if this boundary layer actually oeeu;ies half the tusnel width b/2, i 
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de ` By combining (5) end (6), the ratio ef these tvo pressure losses is 

a a. HM Z 


= E 7) 
a cie: an RI] 
5 path ratio given by (3), ve have 
08 w 
de Geer - pnus z 
aA * AAT us RR 122, 
and VX = 2 na before, then vhen the rette (/ 19 0.056, the ratio of 
frictions] lees to shook lose is 0.628. ihe 


frietional loss and 
-the shock loss is of the seme order of magnitude, 
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These element..y caleulations(mekes it clear thet for the design of the 


nonale ani test section for « superwerodynamics wind tunnel, it is no longer 
possible to separate the compressibility effects and the viscous effect. In 
fact, the concept of boundary layer is also of doubtful v.lue due to the 
extremely small Reynolds amber encountered. Therefore to cctuslly design 
such a nozzle to obtain the mearost approximation to the idesl uniform flow 
it vill be necessary to use the exact Mevier-Stokes e tions instord of 


* D 
[000 g = (ae Ze oy cle. lera = Ae xto an 


tne ip roaizate boundery lsyer eau tions. 
for su«r'erodjnamics, the tavier-Stokes e wii us for no nore exact end 
edd. tlonal corrections aust be +dded (Ref. 2). However, recent !nvesti.a- 


Of course, It miy be argued that 


tions by R. ichemberg (Ref. 3) have shown that these sddition:l corrections 

ere smell in cse of slip-flows concerned here and will not essentially 

alter thé flow jettern. Hence for a first approximation Just like the non- E 
viscous isentro;ic flow as a first approximation for ordinary supersonio 

H Uoisles, we can use the Navier-Stokes equations. The simplest case to be 

cousidered is certainly the axially symmetrie nosrles. If x is tne coordinte 

in the axial direction, 7 the coorainste in the radial direction end u 


und v are the corresponding velocity componente, (| the equations arer 
N E + 426m) ee (9) 
d WA (20) 
N pap = -2% + Irad (T), (u) 
f 


PEL E eer) -8-JuGrdlrh to fttt, fr BH) f D un 


vhere Beni; 
^C t 
Z = dessipition fonction 

T = siresses tensor 


(9), (10), (11) und (12) togetner with the equation of states 
* =R? (13) 


then determines the five unenows u, v, É, f, and T. Of course, the actual 
process of making this celoulation will be extremely tedious and some aj rox- 
imation method of solution may have to be developed. One possibility would 
be to adopt the Karman-Pqlhausem method for boundary layer to this o. ren we 
inthgrate the differential e wationa ence with respect to r and thus only 

try to satiefy the equations fon the Average“ over the cross-section of the 

+ norsle. Tbe "distribution" of u, v over the cross-section will then be set 

in the form of ( ;olynomiesl in r. Initial study im this procedure is already 


mide by 5. A. Seat (Ref. 4) at the mupcestlon of the author. 
For oriioery supersonic diffuser, high efficiency of pressure rcoovery 
Can be generally .chieved by using a long diffuser. Howaver, for supernero- 
dyn-mics wind tunnel, due to the extremely lirge loss through friction, long 
diffusers ure undesirable. In fect, the pressure lose con be reduced by 
using » shortest ;ossible diffuser. 


2. 


Flow Meesuresent 
The .unntities which determine the fiov field ors te out of the 
three virisblesf ,9 , T and the velocity coz,orents. The „usntivtesh, P , 
Tre related by tne equstion of states sad therefore only two i» necessury 
for the detersination of all three. Generaliy for wind tunnel work, the 
quantities actully measured are f, fand f, the magaitude of the velocity. 
or the urement of pressure, „ afanomet r is uced. For ordinary 
pressure, ons uses a fluid = ester filled witn water, zlconol or meroury. 


However for tne extremely low pressure encountered in the suparnerodycante ‘ 
flow, sone other form of manoueter is necessary. One of the aost weed type 


is the Pirani gage. The conventional fora of Pirani ysuges hau ` pressure 
ensitivity of about 10°? mieron.® It utilises the ch ade of temperaturo 

of u wire ber ted with constant energy caused by a ele pen in the prossure of 
the gus surrouncing it. The temperature chinge 1s measured by the ch nge in 
the resistance of the vire. The wire is located in » smali ch aber tileh 

is connected to the „oint oi measurement by s bole, flusn with the gas stream 
Af static pressure ic to be measured. The question of best design of the 
connecting tube for wich response is studied by 5. 4. Schaft. (Ref. 5). 

To measure the density? » the conventions! trod utilises the differ- 
ence in the velocity of light rays in mediums of different 4anatty. ith 
different optienl errergements, ve hive the shsdowgr&ph method, tio schlieren 
method na tie interferoseter method. However, if the density of the Medium i 
is very low se the case of superserodynamic flow, the sensitivity of thee 
methods become extremely poor. For instence in c se of schlieren method, the 
percentage change in illuminction I by passing through - region of thicknese 


D 
Jo leren = 1a a. Hg. one seelen = OO EP miiran 


b ts given by 


I 
F (Fee If fl au 
where Ô is the air density at 32°7 and 1 esl here ,ressure, and 4//2% 
is the density gradient normal to the light rey. fand ^re the local 
length and the corms] unobscured width of the light source imaie erkenate- 
War to the fe edge. K is fester of order 1, determined by the particular 
o;tical „atn used. Therefore the sensitivity of the scilieren xethod de- 
eren:en with the factor ( f/f, ). Some imporvement c-r be mnde by altering 
the quantities f sud & but pruetical limitations and dif!raction difficulties 
do not eilow tie increcse of sensitivity to setisfactory values, 

A new »pproach to this ,roblem density mesnsureaent is the meticd of 
ebsorption. It is found for instance, that oxygen at low pressures shove 
a strong absorption band at vave lengths around 1% j or ultra-violet 
Light. The percentage :bsorption is proportional to tie number of molecules 
that meet the light rey and is, therefore, proportional to the density of 
the gas. The measurement is then similar to thet of toe interferoneter 
method where the density of gas is determined. A similar method is the 
utilis:tion of toe after-flow of nitrogen. These methois re nov being 
studied by R. A. Evans (Ref, 6). 

The conventional method for the meusurement of velocity i» through the 
use of dynamic pressure rise in a "itot-tube. A straight forv:rd application 
of this method is, hovaver, difficult for rurefied geses. The formule used 
is, however, b.sed upon the neglection of viscosity effects, But for rere- 
fied gases the viscosity effect 16 of greet Importence aè ‚olnted out in 
the previous section. Then the dynamic pressure would be „ulte different 
then that given by the usual formula. To entiante this effect, let us 
consider the case oí low Mach number so that compressibility effects cun be 
neglected. Then as « ‘iret approximation, take the flow Field around the 
Pitot-tube as tnat of a source af strength) in non-viscous flow of uniform 
velocity U. (Peg T). The "radius" of the tube & is 
Len 


«nd the stagnrtion oint is located at 


Pet Lis 
2 — 5 
3:27—' measuring device to euper:erodynamic flove, oae is naturclly le^ 0 the 


thought of otner avenues of approach, One poasibilit; is the use of hot- 


Rf | P 
Hu” 2 
The velocity introduced by tne source is then 
— 777 5 
By celcul:tiag the viscous stress from this spyroximete disturbence velocity, 
we b.ve for flow "long the axis 
U g LM yy d pa 
Bit U^ 7 
Hence if d is the stagaetion ‚ressure end f the static cressure, 
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Ct ar Si 
Yor rorefied genen, the valus of X/«[/ or the reciprocel of the Reynold's 
nuaber of tha Pitot-tube could be of the orcer of unity. Then the dynante 
pressure rive L ie not the ususi velues Ef" vu * value much less — 
than that. V 
ke calculette the-welositx//. 

hen the velocity of flow is high, ve have the added csg lief tien due 
to the shock. Tne conventional Rayleigh formis for Pitot tubes in sujer- 
sonic :lov is based upon the assumption of very thin shock weve chi of 
the Pitot tube. Now the thickness of the shock ia proportionti to the mean 
free jath. Hence ia rurefied Deep, the thickness of the shock vill be so 
Ancreasec ve to cause interference vith flow in the neighborhoo! of the 
Pitot-tuse. Tris together with the viscous effect mentioned in previous 
paragraph definitely show the imapplicability of the Reyleigh romul: for 
supersonic velocity of rarefied gases. 
tei the great cos licutions In applying the conventions) velocity 
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J. astez of Flow 
The tvo parameters thet ure directly connected with the flow field 
are the Aeynold's arar Re, defined as 


ie = 


77 
where y'is the xinemetic viveosity,L the typicel lineer disension of the 
bodyj and tne Mach number jf of the free stroan. This is true even for 

slip flows and free solecule flows due to the fact thet the ratio of meon 
free path to the typical dimension can be also expressed in terms of the 
Reynold’. number «nd the Mach ater. 

However, as the pressure or density is reduced, the solid boundery of 
the flow enters ectively into the flow conditions by re,uiring not only that 
‘the microscopic stream velocity be tengential to the surface but that the 
interaction of the molecules «nd the wall be cou red end that the redia- 
tion of energy to and from the wall be taken into account. The interaction 
of the molecules with the vall is so far ex-ressed through the fraction 5 
of molecules that ore diffusely re-anitted f the necowsadetion 
coefficient . (We hive reasons to believe vare functions 
of the tes erature of the vali 7 distribution, 
Therefore the interaction of tne molecules with the wall is the same only 
if the wall teuperecure, the tem, ers ture and the Mech number of the 1 
vun above the wall is the sme. These considerations sees to Indicnte then 
that ‘or the model test to be similar to the prototype, the model pust be 
made of saxc surface moterial ss the prototy;e, and the fluid aust be the 
Baus, eno surtuermore tne foliowlag parameters must be the saner 

(i) Reynolds number Re 
(2) Much number M° f 
(3) free strona tespersture, 7° 


n 


The ridivtion hect loss from the surfece is ecual to TL per unit 
aree. However, if the model is surrounded by the wells of the test ale ber, 
there 1s lee az he«t input due to red! tlon from walls of the test section 
to the model, Lot us call this ,uantity 7, . Then tha net heit loss er 
unit «rer of the eurfsce of the model 10 f- f, . This quantity can be 
rendered ron-dinensional by dividing it by p Ki 7°) - Call tis new 
pıraneter Am , then 
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For the prototy;e, the heat from the walls of the Lost chamber is ubsent 
t but theru may be eoler rediation and tbe radiation from the earth and sur- 
rounding +taos here. (Ref, 9). Denoting tnis mount d : then the 
parameter / for the prototype is 


Pre (36) = 
A= "euer 
In order for the flow to be same also with res.ect to the rhale Uon heat 
transfer, m 


^*^. (Ga) 


Because of the previous conditions on the Reynolds number «nd free stream 
Veaperiture, (27) is the eame as 
"Test. , 
fett e» 
where L, 1a the rie) limecr dimension of the model and L te the typical 
Linear dimension of the prototype. This mecns thet the wall tenperatura of 
j the test chamber must be so controlled that f, s:tisfies (#4). 
f This set of rsther strict similarly rules for model testing in surer- 
serodymenic flows is certainly difficult to satisfy. In what way tho rules 
can be relaxed is the problem of future rese:roh. 


THE “LIMITING LINE" IN MIXED SUBSONIC AND SUPER: 


Hsveeshen Tsien 
California Institute of Teehnology 


It is welleknem that the vorticity for any fluid element is constant 
if the fluid is non-viseous and the change of state of the fluid is isentropic. 
When a solid body is placed in a uniform stream, the flow far ahead of the 
body is irrotational, Then if the flow is further assumed to be isentropic, 
the vorticity will be zro over the whole field of flow. In other vors, 
the flow is irrotational. For auch flow over a 


lid body, it is shown 
by 7. Theodorsen (Ref, 1) that the solid body experionees no resistance. If 
the fluid has a small viscosity, its effect will be limited in the boundary 
layer over the solid body and the body will have a drag due to the skin frio" 
tion, This type of 


entially isentropic irrotational flow is generally 
observed for a stream-lined body placed in a uniform strem, if the velocity 
of the stream is kept below the so-called “critical speed", 

At the "oritienl speed" or rather at a certain mine of the ratio of 
the velooity the undisturbed flow and the corresponding velocity of sound, 


shock waves appear, This phenomenon is called the “compressibility burble", 


Along & shock wave, the change of state of the fluid is no more isentropio, 
although still adiabatic, This results in an inorease in entropy of the fluid 
and generally introduces vorticity in an originally irrotatioral flow. The 
increase in entropy of the fluid is, of course, the consequence of changing 
part of the moohanieal energy into heat energy. In other words, the part of 


fluid affected by the shook wave has much 1: mechanical energy. ‘Therefore, 


with the appearance of shock waves, the wake of the stream line body is very 


much widened, and the drag increases drastically, Furthermore, the acconpnied 


change in the pressure distribution over the body charges the aurodynanie 


moment acting on it and ir oase of an airfoil decreases the lift force. 
A11 these consequences of the breakdown of isentropic irrotational flow 


are ro;orally undesirable in ap lied ecrodynamios, Its ocourrence should be 
delayed ds much as vossible by modifying the shape or contour of the body. 
However, such endeavor wil! be very much facilitated if the cause r the 


ort ter in for the broakdown oan be found first, 


Criterion for the Break: 


G. I, Taylor and C. F. Sharman (Ref, 2) oaloulated the successive 
*pproximations to the flow around an airfoil by means of an electrolyte tank, " 
They found that whai the maximum velocity in the flow reaches the local vel- 
ooity of sound, tho convergence of the successive steps seems to break down. 
| This fact led to the identification of oritioal speed or critical ‘ach number 
j with the “ach number of the undisturbed flow for which the local velooity at 

some point reaches the local velocity of sound, However, there is no the- 

matioal proof for the coincidence of the critical Mach number so defined and 


the breakdown of isentropie irrotational flow, Furthermore, such a definition 


for oritioal "ch number implies that a transition from a velooity less than 
that of sound, or subsonic velocity, to a velocity greater than that of sound, 
or supersonic velocity, does not occur fr: isentropic irrotetionel flow. Un 
4 the other hand, Taylor (Gef. 3) and others found solutions for which such a 
transition occurs. Furthermore, A. . Binnie and S. O. Hooker Let, 4) have [ 
shown that at least for the oase of spiral flow the method of succe sive 
approximation ie a convergent one even for supersonic velocities, With these 


facts in mind, one can conclude that the idertifioation of critical speed 


with local supersonic velocity cannot be correct. 
Taylor's investigation on the spiral flow (Ref. 3) indiontes that there 
is a line in the flow field where the maximm velocity is reached and beyond 


which the flow oannot continue, W, Tollmien in a subsequent paper (Hef, 5) 


called such lines as the “limiting lines". The velocity at the limiting 


line le never subsonic, However, the true characteristics of such limiting 
lines and their signifioaroe were not investigated by Tollmien at that time. 
Recently, F. Ringleb (Nef. 6) obtained another partioular solution of isen- 
tropic irrotatiorel flow in which the maximum velooity reached is approximately 
twice the local sound velocity. For this flow also, a limiting line appmred 
beyond which the flow cannot continue. Furthermore, he found the singular 
charactor of the limiting line, such as the infinite acceleration and Infinite 
pressure gradient, Th. von Karman (Ref. 7) demonstrated this fact for the 
gere wel two-dimensional flow, He also suggested that the limiting line ia 

the envelopo of the “ach waves (Fig. 1) and thus oan only our in supersonic 
region. ie also took its aooearunoe as the criterion for breakdown of isen- v 
tropio irrotational flow, This rencral two-dimensional theory ma established 

later by both lingleb (er. 8) and Tollmien (lef. 9). Tollmien oorreoted 


some mistakes in Ringleb's paper and in addition, showed that the flow de= 


finitely cannot continue beyond the limiting line, The later fact introduced 
a "forbidden region” in the flow bounded by the limiting line, This physical 
absurdity oan only be avoided by relaxing the condition of irrotationality, 
But am stated previously, for non-viseous fluids, the transition from a flow 
without vorticity to that with vortieity oan only be accomplished by shook 


waves, which at the sane time also cause an inoræse in the ontropy. 


However, before one can conclude that the appearance of limiting line, 
or the «mvelopcof Mach waves, is the general oondition for breakdown of 
isentropic irrotational flow, one must prove that the singuler behavior of E 
limiting lines cre penoral end not Limited te two-dinensioel flow. This is 
the purpose of the prosent paper, First the property of limiting line in 
axially symmetric flow will be investigated in detail. Then the ger aral 
throe-Aimensimal problem will be sketched, These investi” tions confir 
the r-sults of Ringleb, von Karmın and Tollen for mene nore general oases, 
Therefore, ty considering only the steady flow of non-visoous fluids, 
the criterion for breakdown of isentropie Irrotational flow is the appearance 
of limiting line, However, for the actual motion of & solid body, the flow 
is neither steady ror non-visoous, Small disturbances almys our and almost 
all real fluids have appreciable viscosity. The sml] Jisturbanoes in the 
flow introduce the question of stability, In other words, the solution found 
for isentropie irrotational flow my be unstable even before the appearance 
of limiting line, and tends to transform itself to a rotatioml flow involving 
shook waves at the slightest disturbance, If this is the case, the criterion 
concerns not the limiting line, but the stability limit. This problem has yet 
to be solved. 
The effect of viscosity will be limited to the boundary layer if the 
pressure along the surface in the flow direction rever inoroaees too rapidly, 
Then outside the boundary layer the flow is isertropio and irrotationul, 
If the gradient of pressure le too large, the boundary layer will separate 
from the surface. However, at low velocities su h separation only widens the 
wake of the body and changes the pressure distribution over the body. But 
if the boundary layer separates at & point where the velocity outside the 


-5- 


boundary layer ir supersonic, a'ditional effects may appear. The flow outside 
the boundary layer in this oase can be regarded approximately as that of a 


solid hody not of original emtour but of a new contour including the "deed 


water" Carl mm orantad by the separation, It is then icmedintely olear that 


the ideal isentronic irrotatiaral flow around this new contour my have a 
limiting Line. 


lenoe, the actual flow then must involve shook waves. In 


other words, the senaration of boundary layer in supersonic rerim may induoe 


a shook wave and thus extend its influence far beyond the rerion of separa= 


tion, Furthermore, the eteeo adverse pressure gradient across a shook wave 
may wocent the separation, This internation between the separation and the 
shook wnve 15 frequently observed in experimenta, 

The above considerations indionte the possibility of the breakdom of 
Asentropic irrotational flow outside the boundary layor even before the 
appearance of limitine line, Therefore, the lach number of the undisturbed 
flow «t which the limiting line appears my be called an the "upper oritioal 
Mach number," On the other hand, sinee shock uve oan orly ossur in supere 
sonio flow, the Mach number of the undisturbed flow at which loo] zl sei 
reaches tho velocity of sound may be called ns the "lower orition! Mach number", 
The aotun] oritioal Mach number for the nppenranoe of shock wavos and the 
compressibility burble mst lie between these two limits, Ay mrefully 
Aesiming the contour of the body to avoid the crowding together of Wach 
waves to form an envelope and to eliminate adverse pressure rredient nlonr 


the surfoe of the body, the compressibility burble can be delayed. 


Axia! m 
The solution of the exact differential equations for an axially 
symmetric isentropic irrotational flow was first given by F. Frankle (Hof. 10) 


The method is developed independently by C. Ferrari (er. 11). Their method 
applies partioularly to the oase of supersouio flow over a body of revolu- 


tion with pointed nose, In this oase, the flow ut the nose can be approxi= 


F mated by tho well-known solution for a cone, From this solution, the differ- 
ential equntian is solved step by step using the net of ohermoteristios which 

| " are real for supersanie veleeities, In the following inv-stigntion, the 
chief eonearn is not the solution of the partiel differential equation but 
rather the oocurrence and the properties of the limiting line in an isentropie 
irrotational flow, The general plan of attack is that of Tollmein (et, 9). 
However, hore the onleulatian is base! on the Legendre transformation of 
welooit potential instead of the stream funotion. 

If 7 is the magnitude of the veloeity, 2 the corresponding veloeity of 
sound assuming isentrop's process, $ the pressure and f the density of 
fluid, the Bernoulli equatian gives 


f- EE igj” a) 


25 = SE = (mE) (2) 
$- 2072 lit n E 2 Á i 
} In those equations, the subseript o denotes quantities corresponding te 7-0 


and / is the ratio of speoifio heats of the fluid. Let the axis of symetry 
{ de the x = axis, the distance normal to x - axis be denoted by y, and the 

velocity components along these two directions be denoted by u and v rospoo= 
H tively (Fig. 2). The x = y plane is, therefore, a meridian plane, Then the 
| kinematioal relations of the flow are given by the vorticity equation 

Lap 00 Ei 

‘Tire Gout thls paper, partial derivatives are denoted by subscripts. Thus 
EE 4 iy. 


and the oontimuicy equation 

dief E Ae 

aO DE (8) 
Eqa, (1) (2) (3) (4) and (5) together with the relation 2 wer" reif 
the flow comletely. 


To sinplify the problem, a velocity »otentíal / Aefine a» follows is 
Antroduoed ! 


u= %, vez % ei 
Then Pq. (4) is lden ny satisfied and Eq, (5) together with Ede. (1) 


and (2) rives the equation for . 


(- Ey, - 2 4 hy +l- Ey + Er 


The characteristios of this differential equation, to be called the cr 


Wl 


aoteristios in the physical plane is given by flug, » where žag) 
is determined by the following equation 


EA LEAN tli- EI e e 


It oan be easily seen from thie equation that 2 is ral only wen 222. 
Theroforo, the oharaoteristios are real only in supersonic rogions of the 
flow. 

The moaning of oharecteristios in the plysicel plane is irnotintely 
aleur if one omlouintes the relation between the slone of a ohoruoteristio 
and the slope of a stroam line i: the meridian or xy plane, Hy the do^inition 
of the function g(x,y), the value of g ia sero, or constant, alos; a ohar- 
acteristic, Therefore, by writing & quantity evaluated at & certain constant 
value of e parameter with that paranct-r as a subscript, the slope of 
characteristic in physical plane is 

(9) 


ba e 


Along a stronn 1 be, the strom furction // defined by following equations 
do constants 
ET uA 
% gu b-p” ao) l 
Therefore, the slope of a stream line is con 


h- t > 


Das. (8), (9) and (11) give 


d 5 , F tang f- {itl mal © 


whore J ia the Mach anle given by = Jn KE? + Therefore, Eq. (12) shows 
d that the oharacteristios in physioal plene are inolined to the strom lines 
by an angle equal to the uch angle. Such lines are the wave fronts of in- 
finitesimal disturbances and are oalled ach waves, In ther words, ohare 
aoteristios in physical planes are the Mach waves in that plane, There 
are two families of Mach waves inclined syrmetrioally with respect to each 
stream lino. 
If to each pair of values of u and v, there is one pair of values of 
x, y, then x and y oan be oonsidered aa functions of u, v. In other words, 
j instead of taking x and y as independent variables, u, v oan be used as in- 
| dependent variable, The plane with u, and v as coordinates ia called the 
"hodograph plane." An equation in hodograph plane oorresp nding to Eq. (7) 
oan be obtained by means of Legendre'a transformation, By writing 


A fu -f Qs) 
eg, Er oo 


T Then Eq. (7) oan be written as 


0 Eur A. -L EI 3 lod le an 


it is soen that 


The charast “istios of "q. (15) are given by Hlap)-owere f Se the lm 
tion o^ following “fforential equation 


MR A t lE- E LE L 
UI + EEE 


Bq. (16) shows that the characteristios in hodograph plane depends upon the 
Ns values of the derivatives of |, which must be obtained from iq. (15). In 


(16) 


other words, the charecteristies in hodograph plane chan e with the flow 
and are mt a constant set of ourves as those ir two-dimensional problems, 
To obtain the relation between the charactoristios ir physical plane 
and moe in hoðogmnh plane, it is notioed that Fa. (9) can be rewritten 


d de 4 P a7) 
Then Fq. (8) is equivalent to 
(EID, eU) Us, +l- E) -v a») 
However, i: goneral, d. (14) gives the followi:g relation = twee the 


difforantinls of x und y and those of u and vt 


dk o dat Ali, de as) 
4 = i, Mat Ii, d. 
By mere of dete relations, d- (18) oan be transformed into an equation for 
(da), aià (dr), Das transformed equation aan be einplifio! by using 7 


Eq. Za, the final relation is 


(Cer E) k. fef CE ZA hg | 
T J-E) dis fim] =0 


Therefore, if tho first factor of Pq. (20) is not zero, the variations (de), 
and (dr) along a charasteristie in physical plane must satisfy the NEA 


l-g 224 Ki E haufig (tehy vd äis A, 2 


Ki 
8 ie the sane ER for the variations Au, V enar) along a charnoter= 


(21) 


isties in hodograph plane as oan "e seen from u (18) and the folloring ree 
lation obtained fron the definition af 
(Wy: Wah = - f: b (22) 
The transformed characteristios of physical plane and tho characteristics of 
hodogruph plane themselves satisfy then the sane first order difførential 


equation, Therefore, these two types of curves are the sane, In other words, 
the eharaeteristies of hedograph plane are the reprssentatian of Mach waves 


in u, v plane, 
The Limiting Line 
M Eq. (20) shows that if 
E E (23) 


then the transformed differential din for charasteristios of physical 
plano, or "oh wevos is tried. Therefore, if thore ise line in the 
hodograch plane, aleng which the values of the derivatives of / are auch 
that Bq, (23) is true, then this line when transferred to the physical plane 
| will have ite slope equal to that of one family of Mach waves. Such lines are 
called as limiting hodograph in ur plane and limiting line in physical plane, 
Sinos Mach waves only oocur in the supersonic regions, it is then evident 
i that the limiting line must appear in these regions, The signífioenee of 
Í the adjective "limiting" will be made clear as other proporties cf sueh lines 
are investigated. 


Tow tha „.ostion arisosı Can the limiting hodogreph be the aharaateris- 
tio in u-v plane? Along a limiting hodorranh, Sq. (23) rives 


(de) = MP iy = ther PEL Gu) 
7 e E = AER t 2 


where the oer ant / denotes the value alor a limiting hodogranh, Now the 


general difforential equation for / , Bq. (15), is true for the whole u-v 
Dlare, therefore, the equation is still true by differentiating it with 
respeot to u and v, “he results oan be simplified bv using A. (15) itself 
und Eq. (23). Then at the limiting hodography 


IHE) EA e I E En 
= Hed, - Eh t GDE ka 

IC 
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Eqa, (24), (25a) and (25b) are the only available equations involving no 
higher derivative than the third, Or the other hand, the slopo of a char~ 
asteristic in hodograoh plane can do oalouleted by d. (22), 


(EI re 4 (26) 


This equation torother with :q, (16) gives 


IRI ger f] n f EE ef in 


Therefore, if the limiting hodograph is a characteristic, then Kä must 


satisfy E4. (27). However, a simple caloulation shows that it is not even 


possible to obtain a relation between 22 and other quantities not ine 
volving the third order derivatives of | . Henoe, e) does not satisfy 


Eq. (27), In other words, the limiting hodogra h is not a characteristic. 


Transferred to physical plane, this means thet the limiting line is not a 
Mach wave, But as show: in previous paragraphs, the limiting line is every 
where tangsut to one family of Mach waves, Consequently, the limiting line 


must be the orvelope of a family of Ach wives. This property of limiting 
line oun be taker. és its physical def At. 


Linitlag Locograph end the Jtream Lines 


At the limiting hodograph both de. (15) and (23) hold, by oliminating 
onc of the second order derivatives, say Lë » the following relation is 
obtained 3 
= =f 
acce E (20) 
(Ave), TE ike) 


The sign before the redionl in d. (28) oan be either positive or negative 


but not both, hie relation will be used presently to show that the stream 
lines and one family of characteristics are tangent ir uev plane, 
From q. (10), the differential of stream funotion onn be onloulatod 


ERT (e) 


In this equation, y can be replaced by x according to d- (14) aw the 


differontinlo ji end replaced ly the Ailfercutialdy and A. coocrül:g to 


iq. (10). Then 


Ada Aug [U rhat thy u. (onus ved. EL 


Aë 


Along a stream line, PE therefore the slope of stream line 1 hodograph 


f plane is given by 
j gi. her S 
po cru ur 
At the limiting hodograph, Eq. (23) holds, therefore, Fq. (31) together with 
Eq. (28) rives 


x NE > 

(de) x rs (de) 5 
ye Lë A L n -/ 

whore tho sign before the radical oan be either ncyative or positive corres= 


ponding to the sim in Eo, (B6). | 


(s2) 


On the other hand, the slope of the characteristics in hodorraoh plane 
ie determined by “o, (27), By solving fer Z) and simplifying the result 
with aid of "e, (15), 


w £ 
Mi) - Alq e 


as (Ar Zi 


The sign before the radical is either positive or negative corrosponding to 


(53) 


the two fanilies of charaoteristios, By using the oositive sign in conjuno- 
tion with nositive sign in . (20), and similarly for the negative sign, 
2 
IL, 
Zi 0 S86 $8. 
Zt at A £ ba d 
Eqs. (32) and (34) show that the stream lines and one family of oh.raoteriae 


(34) 


tios are tangent to cach other at the limiting hodopraph. ‘hie result is 
the same as that obtained for R1 flow, (Ver. 7, o, 9), “hose 


equations when campared with A, (12) for the elope of ach wi 


in een 


plane yields the interesting result that the streum lines and onc funily of 


characteristios at limiting hodograph are perpendicular to the corresponding 
Mech waves at tho limiting line, 


GI - - A (38) 


Eq. (32) gives the following equation which holds at the limiting hodograph 


Eli, + seen l-E m 


CC 8 


Sinee 


coordinates in u, v planer 
4-7 C, v7 Zog 
where Cis the mele between the velocity vestor and Maxie, Then "a. (38) 


takes the form 
( y, M . Zh ILAN =0 (37) e 
This oan bo romrded as the equivalent to iq. (25) for def ning the limit 
hodograph. Similar relation existe for two dimensionel flow, (er. 7, 8, 9) 
L j Along a stream line, tho ratio between Hr}, ana 22 do given by 
Eq. (31). Sy eubstituting this ratio into 4. (19), the E ay 
id), alog a stream line is given as 


n da] 


d 7 , -v dy l Air, be 


At the limiting line, . (23) is satisfied, Then tq. (38) shows that at 
the Limiting line, the stream line has a singularity. Or, rore plainly, Ke 


and KO at these points are infinitesimals of higher order than Ke 
and dr N writing S for the distanoe measured along a stream line, 
Eq. (38) gives immediately 


ES Vit toe 
era H 


éi m 


Therefore, at the limiting line, tho aocelerntions along a ren line is 
Aröinitely Inge. furthermore, sinse the prossure grudiont C A alorg 
a a*5rem line is 

hh -91h =-pfulah tv ly] (a) 
the pressuro gradient al the liniting line is also i finitely large. 


Such infinite acceleration and pressure gruiie;t loud o.c Lo suspect 


TT 


‘that the fluid ie thrown back at the liniting line. In other words, the stream 
lines are doubled buok at this lino of singularity, To L;vesti ate mother 


this te true, tho chamater of the relation yy dry — tap =O donna 
stream lino hes t» be determined. If the dorimtive of thie exrension along 
^ strom Line Am not rero, then d Par Ae. has only a simple zero 
at the intersection of the liniting line and the strem Line, Tonsementirs 
the Alffarantiola (diy and Kä will chango sim by raire throush the 
Linttine hod meh in wer plane along a stronm Mire, Hanos, the straan lines 
WAIL edle baok and form e ouso at the limiting line, The derimtire of 
Xd "c Along, the etrean Line mar be zz: ted with the 


aid of fq. (30) 


LIÉ GLA va Tor 4] = che eub * A ler 


ech d 2 d 
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The expression on the right of Eq, (42) cannot be reduced to sero by the 


(a2 


available relations which consists of Eq. (25), Sa, (16) and differentiated 
forms of Fd. (15). Therefore, the expreasion conserned generally only has 
a simple s ro nt the limiting hodogranh and the str an lires are doubled 
baok at the limiting line, It will be shorn later thnt there ie no solution 


nossible hey the limitinr Lire, Hanes: the nane limiting line, 


Keel, of Characteristios i: lod len: 


Sinee the limiting line is the envelope of the iach wav 8 in physical 
plano, it is interesting to see wicther there is also ezvolope for the 
oharasteristios in hodograph plane. The chareoteristios in u-v plane are 


determined by lq. (26). The envelope to then cen bo found by oliminating 


4 4 betw:en d. (26) and the following equation 


DN u) _ fue Y =0 (43) 
IRL E xii) TE as 
which is obtained by equntire to tero the partiel derimtire of Fy. (26) with 
roapect to &) . The result oan be simplified by “q. (15), and then it ie 
sinply 

j= uto? „ 2 u 


a " E (44) 
This is watisficd by either 


8. 
PL T ai (48) 


The first coriition, Eq. (45), when substituted into Eq, (26) gives 


GL së (n 


which shows that the oirole of maximum velocity oorresumding toe = 0, is 


cond omdie 


the envelope to the characteristics in hodograph plene. The 
tion, iq. (48) is tho spurious s«lutioa, sinee generally the chamoteris*io 
at J= a doos not tangent to the sirole 7 = a, Henoe a = 0 is the only 
envelope. 

‘The Lines of constart velooity in hodograph plane ars cimly oirolos. 
Therefore de. J TE 
(da a x (48) 
By wenne of this rolation wid d. (19), the slope of tho lines of ooistant 


velocity is given ns 


EA 2 v Au, Aer. (49) 
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This equation together with Eq. (30) gives the following interesting relation 
ol. - E 7" 60 
2 75 (hy “= 
In other words, a line of constant velocity in physioal plane is perpendi- 
oular to the stream line in hodograph plane at corresponding point. 


The Lost Solution 
Throughout the provious oaloulation, the possibility of using the 
Legendre transformation is assumed, This requires cut for omah pr of 
waluos of u, v there is one and only ono pair of values of x, y.  louovor, 
it is not niways truo, it is possible to have a number of pointy i. Zio 
physical plane nüvi - the suae vnlue of u and v. If chis is tio cuss, thea 
evido:tly it is imosasivle to solve for x and y from the gadr of ‘motions 


u= u (xy), v * (x,y). Xatnemtioally, the situation is expressed by 


saying that the Jaoobin /d, v)/als, 7] /) vaniches in the physical plene, Or 
, 4 = (51) 

However, this is also tha condition for a funotional relation between u and 

Ye . b. „ Y en be expressed as a function of u, In other words, u and v 

i ere not independent. ienee 1f a solution is "lost" or not inoluded in the 

family of solutions allowing Legendre transformtion, then for that solutia, 

y= víu) (52) 
It ie sen that “q. (51) is then identically satisfied, 
By eliminating P from the Seite equation, we obtain 
DEIER E co (83) 
This equation oma be rewritten in the following form by using Eq. (52) 


J- 4. & [« K 2 A P . „ 09 


The vortioity equation, Eq (4) ean be expressed as 


ZA - Myo (65) 


From Zon, (54 ) and (55), one oan solve for % and 4 . The result 
is P à eer 

KREE HCH E x (56) 
Lt EE El] 9 = EE am 


By differentiating the first of . (B8) with respect to y, the second with 
respect to x, the following relation is obtained by subtraction! 


D 2 DA +7 4 = (57) 
Therefore 
T or / E A Ea 


where Hy) an undetermined function of y, However, d. (86) anows that 
for lines of constant values of u where Ji = 4 (de), + 4 (45), =. 


Ec d = Constant (59) 


Q 


Henoe, lines of constant values of u and v are straight lines, ‘This re- 


striction reduoes the function fig) sa Eq. (58) toa numerical constant, Put 
Anne ta. (58) an then 
n Kr (60) 
7 p 


Ma 
Therefore lines of constant value of u and v are radial lires passing through 


the point x = K. Thus the lost solution is nothing but the well-knom solution J 
for the flow over a conical surface, 

From Ya. (59), it is seen that lines of constant velocity are perpen= 
dieuler to the tangent of uev ourve at the corresponding points, By subeti- 
tuting the value of + fron Eq. (57) into Eq. (56a), a relation between u 
and v is obtained: 

van , EJ v EI € 

This is the differential equation for the determinirg the hodograph represente N 
ing the flow cvor a cone. ig, 3 shows the hodograph for a cone of 50° - 
vertex ang e and with a velocity at surface of cone equal to 0.35 o. o is the 
mxinum volooity, i 
drawn from data given by 0. I. Taylor and J, #, Macooll (Hef, 12). 

It may well be mentioned here that the lost solution for tho axially 


+» the value or 7 eorrespording to a = 0. Fig. 8 is 


symmetric flow is not limited to supersonio velocity only as is the oase for 
two=dimensiorel flow. In fact, Wyler and Maoooll show that for mall for- 
ward velocity of the cone, supersonic velocities cour only just after the 
head shook wave. The velocity decreases as the surface of the cone is 
Finally, it becomes subsonic for bol te near the surface of the 


approached, 


sone. Fig. 4 shows a few exemples taken from their calculations (ber. 12). 
The dotted curves in the figure are the Mach waves. The dotted straight 
lines are the boundaries between the supersonic and the subsonio rorions. 
Furthermore, spark photographs of conical shell ir actual flight taken by 
Meoc 11 (Verf. 13) do not indicate the presence of shook waves in rorions of 
flow where such transition from sup.rsonic to subsonio velocities is expected, 
Therefore, at least for this particular type of flow, a smooth tmnsition 


through sonio velocity actually takes place. 


Continuation of Solution Beyond the Limiting Lire 
Sinos it is shown ir n previous paragraph that the stream linos are 


fenernlly turned beck at the limiting line, the question arises: Is it 
possible to continue the solution beyond the limiting line? Of course, 
there arc two mys of oo timing the solutions The new solution de joined 
either smoothly to the given solution at the limiting line or with a dise 
continuity, As shown before, the Limiting line is the envelope af one fanily 
of tho Nach waves, thon at every pol t of this line its direotion differs 
from that of stream lino by an angle equal to the “ach anglo. Dut the ‘ach 
angle im ot zoro except at points where the velocity of fluid hos reachod 
the maximun velooity and the ratio Zo . Tnerofore, the limiting line 
generally doos not coincide with the stream line, an! the cinoontinuity at 
the junotion of the solution at the limiting line ot be that of a vortex 
sheet, The only other type of ^iscontimuity is the shock wave. However, 
the angle betwen the limiting line and flow „ret on is equal to “ach 


angle, Then accor*ing to the result of the theory of shook wmvrs, ths e- 


continuity across such a line vanishes, In other words, thore cannot be 


disconti. it at the limiting lire, Therefore, it is impossible to joi a 
new solution at the limiting line with a disomtinuity, 


to the 


oond possibility of joining a new solution moothly at the 
limiting line, it ic seen that the flow beyord the limiting live must be 
irrotnut'onal a % isentropic since the limitinr line cannot be a shook wave, 
There arc only two tpos of isertropie irrotational flow, namely, one that 
allows the Legercre transformation and one that does not, the lost solution. 
Inv stigate the second alternative first. If the solution beyond the limiting 
line belongs to the so-called "lost solution", then since the junction at the 
limiting line must be smooth, the values of u, and v at the limiting line 
must also satisfy the :q. (61). But, the slope EH at limiting line is 
piven by Iq. (24). The second derivative (de) will then involve the 
fourth order de. watives of |. Besides these expressions, the available 


relations aro de. (23), (15), (26a) (25b) and three more equations obtained 


by differentiating "qs. (25) with respect to u and v, However, it is still 
1 possible for . to satisfy an equation like Eq. (81) where no derivative 
of Bi appears, lenoe, the limiting hodograph does not satisfy the equation 
for lost solution. In other words, the "lost solution" cannot be used to 
continue the flow beyond the limiting line. 

The only remaining possibility is to oontirue the flow smoothly by 
another solution obtainable by Legondre transformation, Smooth continuation 


means that th 


valves of u, v and @ must be the sane at the junction, the 
limiting line. Since shock waves do not appear, isentropic relations still 
hold, The density ¢ is determined by velocity only. The value of u, and 
v are determine’ by the coordinates in hodograph plane, The »osition of the 


limiting lie i physical is determined by Yy , + Therefore, the 


problem oan be stated as follows: At a oertain given curve (, ( in the 


hodograph dane, the limiting hodopraph, the values of V, I, are 
given. ‚| is the paraneter along the given curve, It is required to deter- 
mine a new solution of the differential equation E. (15) with these initial 
valuos, First of all, it 1s seen that with the given data, the le^b hand 


sides of tho following equations are givens 
BR) ka E + he . 
ZIL, Aw # + 2 f (62b) 


rororo 7 
Lae dst Flu (om) 
Kell EH) BEL. LRT (om) 


By substituting those values into Eq. (15), the second degree terms reduce 


Ach. LEB) RIS em 


ër 
which is linear in Aw. + Therefore ,[,, an be uniquely determined by 
"o, 18). In other words, with the given data, the second order derivatives 
of | at the riven curve aih), vA) en re determined uniquely, in spite 
of the fuct that the differential equation (15) is of second degree. k. 
Friedrichs and H, Lewy (Ref, 15) have shown that under these circunstences, 
the function | within a region R (Fig. 5) bounded by two charasteristios 
and the given curve de uniquely determined except on additional constant. 
Co.sequently there ca: be only one solution corresponding b the rive datu 
at the Limiting hodograph, However, this solution is the very one which 


gives the reverse flow at limiting line, Therefore, it is impossible to 


eontinue ts solution beyond the limiting line even by Legendre transformation. 
Since all three alternıtivos fail to offer a my of contimuine the 


solution, the limiting line is truly an impossible boundary to cross, In 
other words, the rerion beyond the limiting line is a "forbidden rerion", 
This nhysiosl absurdity oan only be resolved hy the br nkdom of isentropic 
irrotational flow some distance ahead of the liniting line. 


General Three Dimensional Flow 


The methods used in previous seotiors for investigating the axially 
symmetric flow can be easily extended to the general three dimensional 
ease, In the present seotion, this investigation will be aketohed briefly 
and the results indicated, 
Let the three comoonents of veloeity along the three coordinate axes 
x, y and z be denoted by u, v, and w respectively. Then by introducing a I 
velocity potentinl / defined by i 
ueg, v=, vR m 
the differential equation for / of an teentropio irrotational flow can be 
written as (Ref. 7) 


l thet hy Ma rr p, tArm y deal odor fy 


1f for every trind of u, v, v, there is only one triad of x, y, 2, then 


KA 


the Legendre transformtion oar be used, Thus 


Kragen (en 


Ax dg do (m 
The differential equation for » Eq. (66), 18 then transformed into 
a^ [Bt - fir G4 AB - i] ac E) er Al qun) 
+w AB -H*) + dew Ae + fun HE- B6 ) dur (FG-CH) 


where the “ollowing notations are used 
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By analogy with the axially symetrio oase, the limiting hodograph 
surface is defined as the surface in the u, v, w space, or hodograph space, 
where the following relation holds: 
JAH â 
4=|H BF |=0 
la FC 


The prepertion of this limiting hodograph and the corresponding limiting 


(n) 


surface oar be found by considering the behavior of stream lines and charmo- 
teristics at such surfaocs, 


From Fq. (68), the differentials of x, y and £ oan be written an 


dt Adut PH de e & der (n.) 
Y= H du o B de F du em 
(720) 


dt Gdu+F det dw 
Along a stream line, the differentials dx, dy and da must be proportional 


to us v and w respectively. Thus the equation of a stream line in physical 
space is ayy E EM A bey 
4 Tr libe 


where the subsoript // indicates values taken along the stream line, The 


(13) 


equation of a stream line in hodograph space is obtained by eliminating dx, 


dy and ds from d. (78) by Eq. (TR). The result is 
(diy EN EN, (1) 
aurhy+ Jw Autbr+ fw fer Fu +cw 
whore Z 1s the co-factar of A An the determinante A er Eq. (71), Z the 


co-fuotor . 


. Eq. (74) oan be used in turn to eliminate tw of the 


thros diff rentlale du, dv and dw ir the right of Eq. (72). The result is 


D = HAdu 

ER aut fot fw (75) 

dy) — ___vAde 

by = Aut byt fw Wo 
H (150) 

We), J dré £u. 


At the limiting surface, ^| =p as defined by Eq. (T1), thorefore tho stroam 


lines havo a sinmularity there, Similar to the axially symmetric flow, 


the stream lines generally are turned back and form e ousp at this sur Noe. ! 
The acceleration and the pressure gradient are, of oourse, infinitely large i 
at such places. i 


The characteristic surface LE 2) =0 An physionl space is deter- 
mined by the equation 
Clp + fy th JWI RP ew d yh t duns ftir fy (rm 
Sinoe this equation ie a second degree equation, there are two families of 
surfaces passing through each pol. These surfaces are the wave fronts of 
Anfiniteoimal distrubances in the flow and car be e ed the “moh sr fus o. 
The charecteristic surface E l, V, w) = in the hodograph space is deter- 
nined by the equation 


afta c) f letA)4 Mag AF hh, tb heh, HLA] omy 
DENTEN 

e [817 A ale. + 644-F4 Akt] 
t dwu [ FAL + Hf fe - 64 BLA] 

tdw [ GLU FLA" Ha cue] 


By transfo minz . (76) for Mach surfaces to hodograph space, it oan be 
shown that the transformed equation is satisfied ty either the chumoteris- 
tios in hodograph space determined by Eq, (77) or the Liniting hodograph 
deteminol by 2. (Tl), Therefore, here ain tho limiting sur mes 10 he 
envelope of e fanily of “ach sur Nees. 

By using gs. (74) and (77), it im possible to show that the stream 
lines in hodograph space are tangent to the characteristic surfaces at the 
Limiting hodogresh. Furthermore, by using Eqs. (69), (71) and (74), the 
inolinatior of the stream lines at the limiting hodograph oan be oaloulated. 
In mot, uf (%, e (du) de d del, JÈ Uut the following rela 


-£ mei 


This relation is really equivalent to Eq. (32). In other words, at tho limite 
ing hodogmuph, the trolination of the stream lines and characteristics from 
the Z = constant surface i» equal to the "eh angle (Fir. 6). It thus soom 
the break down of rererel steady isentropic irrotational flow of non-visoous 
fluid 1s oonneoted with the appearance of the orvelope of Maoh w vos in 

bs tes! apace and the tangmay of strom lin-s and oharaoteristios in 


hodormah space. 
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Figure Gapticna, 


Limiting Line as the "nvelope of Mach /aves 


Stream Line and Velocity Components in an axially symetrio 
Flow 


Hodorraph of the Flow over a Cone of S09 telf Vertex Angle 
and ^ Surface Velosity 7 equal to 0.366. 


Flom over Čo es of Various Vertex Angle involving Subsonie 
Rerios. — 2 * Half Vertex Angle, 4, = Velooity over the 
Surface of Cone 


orten R where the Solution is uniquely derer? ed by riven 
es at the Limiting Hodograph. 


"kemotrionl zeien eng betwen tren line and Characteristic 
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April 12, 1943 


Dr. Esue-Shen Tsien 

California Institute of Technology 
Pasadena, California 

Dear Dr. Tsien: 


I just received the program of the Meeting 
of the Mathematical Society at Stanford LM E 
on April 24, 1943. I see that you will address s 
meeting on 

c " you have not yet mai 
arrangements concerning the publicstion of this 


manuscript, I should like to suggest that you consider 
the Quarterly of Applied Mathematics as a possible 
means of publication, 


Very sincerely yours, 


Q ge. 


W. Prager 


April 24, 1943 


Dr. Hsne-Shen Tsien 
California Institute of Technology 
Pasaderra, California 


Dear Dr. Tsien, 


I have noted tuat you plan to speak upon "The limiting 
line in mixed subsonic and supersonic flows of compressible 
Jluids" on April 24th. 


Tids sounds very interesting to ne in the light of 
some recent investigations I have undertaken, and I hose you 
will be kind enough to send re a reprint of your lecture when 
At becomes available. 


Sincerely yours, 


May 7, 1945 


Robert B. Dexter 

Institute of "a Aeronautieal Soienoes 
1606 RCA Building West 

30 Rockefeller Plasa 

Now York, Now York 


Dear Mr, Dexter: 


Fnolosed is my paper entitled "The Limiting Line in Mixed 
Subsonic and Supersonio Flow of Compressible Fluids." I am 
submitting it for consideration of publication in the Journal, 


Sinoerely yours 


H. 8. Zeie 


INSTITUTE 
Mts 
AERONAUTICAL SCIENCES i 
OFFICIAL PUBLICATIONS i 
JOURNAL OF THE AERONAUTICAL SCIENCES | 
AERONAUTICAL ENGINEERING REVIEW | 
lees pak indies ger 


30 ROCKEFELLER PLAZA 4 
NEW YORK | 


May 13, 1943 


Dr. H. 8. Tsien 

Daniel Guggenheim Aeronautical Laboratory 
California Institute of Technology 
Pasadena, Calif, 

Dear Dr. Tsien: 

Thank you for submitting your paper "The 
Limiting Line in Mixed Subsonic and Supersonic Flow 
of Compressible Fluids" for possible publication in 
the "Journal", I am putting the paper through the 
usual procedure and will notify you of the comments 
of the Editorial Board as soon as possible. 

Best regards 


Sincerely yours, 


IL vae Q 4 


Editor 


May 17, 1945 
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INSTITUTE 
or THE 
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OFFICIAL PUBLICATIONS 
JOURNAL OF THE AERONAUTICAL SCIENCES 
AERONAUTICAL ENGINEERING REVIEW 
1808 RCA BUILDING west 
30 ROCKEFELLER PLAZA 


NEW YORK é 
July 22, 1943 

u 
Mr. Hsue-shen Tsien 9 
Guggenheim Laboratory 
California Institute of Technology 
Pasadena, Californie P 
Dear Mr. Tsien: 


Enclosed is your paper, "The 'Limiting Eine. in 

Mixed Subsonic and Supersonic Flow of Compressible Fluids", 
which you submitted for publication in the "Journal". This 
paper was returned to us by the War Department with the 

following comment: "Objection is taken to unclassified pub- 
lication of the article entitled, 'The Limiting Line in Mixed 
Subsonic and Supersonic Flow of Compressible Fluids'. This 
paper contains information that would be of aid to the onemy." 


Your paper is an excellent one and was highly 
recommended for early publication by Dr. Dryden. I regret 
that we are unable to publish it at this time. Possibly, if 
you do not find means of publishing it as a classified report, 
you may be willing to submit it at some later date for pub- 
lication in the "Journal". 


Thank you for submitting the paper to us. I hope 
tinue to send future papers for possible publica- 
tion. I should appreciate your acknowledging the return of 
this material to you for our records. 


Best regards. 
Sincerely yours, 
Secretary 

RRD:BK 

Enc. 


September 14, lud. 


Mr. R. G. Robinson 

tional Advisory Committee for Aeronautics 
1800 Mew Bo Ave., Dupont Circle 
Washington 25, D. C. 


Dear Mr, Robinson: 


About a year and a half ago, 4 submitted to NACA a paper entitled 
"The Limiting Line in Mixed Subsonic und Supersonic Flow of Compressible 
Fluids" at the suggestion of Dr. von Karman for consideration of publica- 
tion as a classified report. Six months ago, Dr. von Korman ms in 
Washington and has talked about this matter with your staff, nd found 
that the paper was mislaid. However, Dr, von Karman told me thet it was 
decided to publish this paper at an early date. 


Since another six months have passed, 1 wonder whether you will per- 
mit ne to inquire again about the publication of the paper. âs this paper 
discusses a subject of great current interest, its carly publication should, 
perhaps, be considered. 


1 apologize for troubling you with such a trivial mutter, ond 4 remain 


Very sincerely yours, 


Baue sten Tsien. 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


1500 NEW HAMPSHIRE AVE., DUPONT CIRCLE 


— Wasuneron 25, D.C. 
vou September lj, 19d "mem — f 

— anii 

re 


Dr. Hsue-Shen Tsien 
California Institute of Technology 
Pasadena, California 


Dear Dr, Tsien: 


With your letter of August 18, 1913 you forwarded 
to the Committee a copy of your paper "The 'Limiting 
Line’ in Mixed Subsonic and Supersonic Flow of Compress- 
ible Fluids" for consideration of its release as a 
classified Committes report. 


Study of this report by the Committ 
staff has been completed and I am p 
that your paper will be published as & Restricted Tech- 
nical Note of the Committ: If 1t appears advisable 
to refer the report to you in its final form prior to 
publicstion, and I suspect that this may be the case 
because of the amount of mathematical work included, 
the Committee would like to refer a proof copy to you. 


s technical 
to advise you 


Thank you for your interest in making the paper 
available to the Committ staff and in offering it 
for publication. 


Very truly yours, 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


Hage — > EA 


G, W. Lewis 
Director of 
Aeronautical Research 


HER 


September 16, 1944. 


Dr. à. u. Lewis 
Director of Avronautical Research 

hs Mon) Advisory Committee for Àeronüutios 
1500 New Hampshire Ave,, Dupont Circle 
Washington <5, D. 0. 


Dear Dr. Lewis: 


It is indeed a happy coincidence that only two 
days ao I wrote to Mir. R. G. Robinson injuiring about the pub- 
lication of my paper "The “Limiting Line" in Mixed Subsonic and 
Supersonic Flow of Compressible Fluids" and then I received your 
letter of September 14, 1944 confirming its publiostion as a Re- 
stricted Technical Note, Following your sugzestion, I would like 
to sec the report in its final form as the paper involved con- 
siderable mathematical work. 


I thank you for informing me, and + remain 


Very sincerely yours, 


Har pl Heue-shen Tsien, 


September 23, 19l, — — (i5 


„ Hsue-Shen Tsien 
California Institute of Technology 
Pasaden: California 


bear Dr. Tsien: 


1 now have more complete information on your 
paper "The 'Limlting Line! in mixed Subsonic and 
Supersonic Flow of Compressible Fluids" and in trans- 
mitting 1t T wish to acknowledge your letter of 
September 10 and one of September 1 addressed to 
Mr. Robinson. 


A mimeographed proof copy will be ready approxi- 
mately three weeks from now and will be forwarded for 
your review before proceeding with the duplicating 
process. The fommittee will be pleased to have your 
comments at that time and the T will be published 
and released immediately thereef . 


Very truly yours, 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


Wl er EL. 


For G. W. Lewis 
Director of 
Aeronautical Research 


RESTRICTED REGISTERED M 


un NATIONAL ADVISORY COMMITTEE 


October 17, 1944. 


Velifornia Institute of Technology, 
Pasadena, California, 


Dear Doctor Tsien: 


With reference to my letter of September 23, there is forwarded 
herewith, by registered mail, for final check before release, a proof 
copy of the restricted Technical Note No. 961 entitled "The ‘Limiting 
Line' in Mixed Subsonic and Supersonic Flow of Compressible Fluids," 

of which you are author. 


It is requested that any necessary changes or corrections be 
indicated in red and that desirable changes or corrections be in- 
dicated in black in the enclosed copy of the Technical Note. 


As it is desired to release this Technical Note in tho very 
near future, I would appreciste your returning the corrected copy, 
together with the enclosed original material, ae promptly as pos- 
sible. 


Tery truly youre, 


Aeronautical Research. 


Enc. Proof copy and orig. text and figs. 


RESTRICTEU 


October 26, 1944. 


Dr. d. W. Lewis 

Director of Aeronautical Rescarch 

Nutional Advisory Committee for Aeronautics 
Rs.hingtor, D. C. 


Dear Dr. Lewis: 


Follosing your request in the letter of Oe ober 17, 
1944, 1 am returning the corrected proof end origin:l manu- 
script of my paper entitled "The ‘Limiting Line! in Mixed Subsonio 
end Supersonic Flow of Compressible Fluids". There are only two 
short correctione necessery and they appesr on puge 6 and page 24 
of the proof copy. 


It is indeed en honor to publish my puper «s a techni- 
cal note of the Comittee and I deeply uppreciste it. 


Very truly yours, 


Hsue-shon Tsien. 


